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ABSTRACT 



Liquid crystals, a material that exhibits brilliant changes in 
color over narrow temperature bands, have been successfully used to 
study the temperature field produced by spherical and cylindrical 
radio frequency surgical probes. An agar-water mixture was used to 
simulate tissue. Experimental temperature data were obtained by 
photographing the liquid crystal material through the clear agar. 
Comparison of experimental data with an existing theoretical solution 
was excellent, agreeing within 8% for the spherical probe. However, 
the cylindrical data and cylindrical theory agreed only to within 30%, 
the large discrepancy being attributed to edge effects. 
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NOMENCLATURE 



Symbol 

A 

0 

A 

b 

C 

°C 

D 

E 

f 

6 

6 ' 

1 
I 
i 

k 

L 

M 

m 

n 

R 



Description 

Surface area 

Wavelength 

Number of flow tubes 

Volumetric heat capacity 

Temperature in degrees Centigrade 

Diameter 

Electric potential 
Frequency 

Nondimens ional power 

Corrected nondimens ional power 

Current density 

Radio frequency current 

Coordinate of nodal point along 
distance axis in finite 
difference formulation 

Thermal conductivity 

Length 

Stability criterion for finite 
difference formulation 

Slope of nondimens ional resistivity 
curve 

Coordinate of nodal point along 
time axis in finite difference 
formulation 

Nondimensional radius to an isotherm 
Ten ohm-ten watt resistor 



Typical Units 
cm 2 

Angstrom 
cal /cm 3 °C 
cm 

volts 

hertz 



amps/cm 2 

amps 



W/m °K 
cm 
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r Probe radius 

o 

r Radius 

r. Isotherm radius 

T Temperature 

t Time 

V Voltage 

a Thermal diffusivity 

y Slope of dimensional resistivity 

curve 

n Ratio of speed of light in air 

to that in water 

e Nondimensional temperature 

A Wavelength 

p Resistivity 

Nondimensional resistivity 
x Nondimensional time 

<j> General potential 

oj Uncertainty 

v Del operator 



cm 

cm 

cm 

°C 

seconds 

volts 

cm 2 /sec 



Angstrom 
ohm- cm 



SUBSCRIPTS 

0 Initial conditions 

10 At the ten ohm resistor 

cct In the circuit 

1 Of the ith component 

p From the photographs 

cyl With reference to the cylindrical probe 

sphere With reference to the spherical probe 
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NONDIMENSIONAL QUANTITIES 



e = 




Temperature 




G 

0 



I 2 P r 2 
o 0 

A 2 k T 



Time 



Power 



R = 



r 
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Distance 




1 + me 



Resistivi ty 
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I. INTRODUCTION 



In neurosurgery, radio frequency current has been used to create 
lesions, thus alleviating conditions such as acute pain and Parkinsonism. 
A current emitting electrode, similar to a hypodermic needle, is inserted 
into the brain through the skull. An inactive electrode is placed else- 
where on the body in an area of high conductivity. Radio frequency 
current, at a low power level, is passed between the electrodes and 
electrical energy is transformed to heat in the tissue with the result- 
ant rise in temperature destroying the tissue. 

To accurately use the radio frequency technique, a surgeon must 
know the exact location of the electrode (probe) tip in the brain, the 
volume of tissue to be destroyed, the thermal and electrical properties 
of the tissue [1], and the temperature near the probe tip. Thermocouples 
and thermistors attached to the .probe have been used to monitor probe tip 
temperatures [2, 3, 4]. However, these temperature sensors indicated the 
temperature of the probe, not the tissue temperature itself. Addition- 
ally, the sensors were nearly as large as the probe tip and may have 
distorted the thermal field by acting as thermal sinks. Without accurate 
knowledge of tissue temperature, tissue boiling occurred with attendant 
gas formation and tissue carbonization [4]. These effects resulted in 
irregular shaped lesions and unpredictable damage along the entrance 
path of the probe into the tissue. 

The objectives of this thesis were twofold. First, to experimental ly 
determine the thermal field produced by a radio frequency probe embedded 
in a resistive agar medium. Liquid crystals, a material that indicates 
temperature through a color change, were used as the temperature sensors. 
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The second objective was to compare experimental data with results pre- 
dicted using a theoretical heat transfer model of the radio frequency 
probe [5]. The general theoretical model included the effects of tissue 
metabolism and blood profusion. These two effects have been excluded 
from the present investigation since the thermal field was generated in 
an inert gel-like substance. 

The temperature fields produced by spherical and cylindrical radio 
frequency probes were successfully studied with the liquid crystal 
material. , The crystals, with their low thermal capacity [6], produced 
brilliant, undistorted, continuous displays of several selected isotherms 
surrounding the probes. Agreement between the theoretically predicted 
and experimentally measured values of temperature for the spherical 
probe was excellent (8%). Agreement between theory and experiment for 
the cylindrical probe was not good (30%). The large differences observed 
in the cylindrical case are attributed to edge effects. The theory assumed 
an infinitely long probe while the experiment employed probes with length 
to diameter ratios of 4, 10, and 20. The edge effects were clearly dis- 
played by the liquid crystals. 
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II. HISTORICAL BACKGROUND 



Use of radio frequency current to create physiologic lesions is not 
new. Carpenter and Whittier [7] describe early investigations into the 
uses of radio frequency current as a surgical tool. Aronow [1] was the 
first investigator of recent time to study the problem of radio frequency 
current in tissue. Lesion formation using radio frequency current is 
accomplished by inserting into the tissue a small hypodermic needle which 
is electrically insulated except at the tip. This electrode acts as the 
active, current emitting electrode. A larger, inactive electrode is 
placed elsewhere on the body. The path of current conduction, i.e. the 
body, acts as a leaky dielectric by dissipating electrical energy in the 
form of heat [1]. This Joulean heating effect results in a local rise in 
tissue temperature. When tissue temperature is increased to 55°C, tissue 
destruction occurs. Accurate prediction of when tissue in the probe tip 
region has reached this temperature has been a major problem with the 
radio frequency probe. 

Excessive tip temperature causes gas formation and tissue carboniza- 
tion which results in poor prediction of lesion size and additional 
damage along the path of electrode insertion [4]. Studies with thermo- 
couples and thermistors have been conducted to correlate electrode 
temperature with lesion size [2, 3]. Since heating takes place in the 
tissue and not in the electrode, the sensors inside the probe tip indi- 
cate only that temperature rise caused by conduction from the tissue back 
into the electrode. This temperature rise may be less than the maximum 
tissue temperature. Therefore, reliable prevention of excessive tip 
temperature has been difficult. The same type of tissue temperature 

measurement problem occurs with the direct current probe, [5j. 
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The direct current probe has a more serious liability than the 
tissue temperature measurement problem. Direct current has preferred 
paths in the body which results in irregular lesions. Death may occur 
if the direct current finds a path to the heart. Preferred paths in 
the body do not appear to be a problem with the radio frequency or the 
high resistance probe. 

The high resistance probe consists of a high resistance iron-nickel 
alloy wrapped around a nonconductive rod. Electrical current is passed 
through the wire causing the probe to be heated. Groff's study of the 
high resistance probe [8] demonstrated the feasibility of producing 
controlled, predictable heated regions in agar. In opposition to the 
radio frequency probe, where heating occurs in the tissue, the high 
resistance probe itself generates heat. Thermal energy is then trans- 
ferred to the tissue by conduction. Prediction of the magnitude of the 
temperature field produced by a high resistance probe thus requires 
knowledge of. only the thermal properties of the medium. Another advan- 
tage of the high resistance probe is the fact that a surgeon may check 
for proper probe tip location before forming the lesion. 

Probe tip location may be verified by raising the tissue temperature 
to 40-49°C using low current levels. Increasing temperature into this 
range blocks nerve functioning but does not cause permanent damage. This 
technique, equally applicable to the high resistance and the radio fre- 
quency probe, is especially helpful in treating Parkinsonism. As the 
temperature increases, shaking diminishes and stops when the probe tip 
has been properly located. Lesion formation is then accomplished by 
increasing the local tissue temperature to 55°C, being careful to prevent 
excessive tip temperatures that result in tissue boiling. 
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While the high resistance probe requires knowledge of only the 
thermal properties of tissue, the radio frequency probe requires 
knowledge of both thermal and electrical properties. Geddes and 
Baker [9] discuss difficulties in obtaining electrical properties of 
tissue. Since tissues are composed of cells, different properties 
may be exhibited depending on the cell orientation. Other factors 
influencing the electrical property values include age, history, and 
environment of the tissue. Tissue properties measured IN VITRO 
differ from those measured IN VIVO. Thermal properties, on the other 
hand, are reliably predicted based on the tissue water content [10]. 
The prediction of lesion size based on the magnitude and duration of 
radio frequency current and tissue electrical properties is, there- 
fore, more complex when using the radio frequency probe than when 
using the high resistance probe. 
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III. THEORETICAL DEVELOPMENT 



A. SPHERICAL PROBE 

Consider a spherical probe, of radius r , which continuously emits 
radio frequency current, I, into a medium that is homogeneous and infinite 
in extent. Assume that the medium possesses electrical resistivity, p, 
thermal conductivity, k, and volumetric heat capacity, C. Gengler [5] 
has derived the basic equation which accounts for the thermal effects 
produced in the medium due to the Joulean heating effects of the radio 
frequency current: 



In the present investigation, the experimental test medium, an agar- 
water mixture, was found to exhibit an electrical resistivity that 
varied linearly with temperature: 

P,= P 0 + y(T - T o ) 

where the subscript indicates an initial, reference condition. 
Substitution of Equation (2) into Equation (1) yields: 



Assuming that k and C are constant, Equation (3) may be normalized by 
introducing the following non-dimensional quantities: 




( 1 ) 




(3) 
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o 

Equation (3) yields: 



1 

w 






30 

9t 



(4) 



The boundary and initial conditions are: 

a) at the probe surface, there is no heat conduction back into 
the probe, 




(5) 



b) the temperature field is undistorted at large distances 
from the. probe, 



R , 0 ->■ 0; 



( 6 ) 
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c) the medium is initially at uniform temperature, 

X = 0, 0 = 0. ( / ) 

Equations (4), (5), (6), and (7) may be cast in explicit finite dif- 
ference form as [5] 







(4a) 



mG e? 
+ — — 
( i aR) ^ 




6 



n 

i 



Ax 



e 



n 

i*-l 



i*+l 



(5a) 



N-l 0 N+1 



(6a) 



6 ? = 0.0 



(7a) 



where i and n define the coordinates of a nodal point in an R, x 
network such that the nondimensional temperature at time x = nAx 
and radial location R = iAR is defined as e^ 1 ; i* is at R = 1 ; 
i = i*, i* + 1,...,N-1,N. Solving Equation (4a) for yields: 
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Define M = At/aR 2 , then 




Since an explicit finite difference scheme has been employed, a stability 
criterion must be stated to limit the size of the time step based on the 
value of incremental distance (aR). Discussions of stability may be found 
in References 11 and 12. The stability criterion for this problem becomes 



1 



m -»• 



MmAR 2 G o 

TW 



> 



0 



Solving for M and assuming the worse case, 



~ mAR 2 G n 

d " iw 

Since M was originally defined as At/aR 2 , choosing a value of AR 2 fixes 
the value of At. 

Computer program I (Appendix E) was used to solve Equation (8) subject 
to boundary, initial, and stability conditions (Equations 5a, 6a, 7a, and 
9). Reference 13 states that the error in an explicit difference approxi- 
mation is on the order of aR 2 . By varying the size of the aR increments 
and plotting temperature versus aR 2 , it was possible to extrapolate the 
exact value of temperature. Figure 1 represents this technique, known as 
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Richardson's technique, for the spherical probe. It can be seen 
that using a value of &R = 0.1 yielded theoretical results within 
2% of the exact value. 



B. CYLINDRICAL PROBE 

Paralleling the development of the basic equations for the spheri- 
cal probe, the equation for the cylindrical probe in nondimensional 
form is: 



1 3 ( R l) + G 0 

R 3R + R 2 ’ 37 

72- y. 2 

where G = p o o 

o 

kT 

0 

1=1= I 
A 2irr L 



( 10 ) 



Substituting G = G q (1 + me) into Equation (10) and transforming to 
finite difference form yields: 

I mG M \ G M 

+ ejll - 2 M + + -fr (11) 

The appropriate boundary conditions are: 
e n - A n 

i*-l " i*+l 



ef 1 = M 



1 - 



2i 



e i-l + + 2T I 0 i+l 



e n - e n 

0 N-1 “ 0 N+1 

6 °. = 0.0 
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Figure 1 

Richardson's Technique 
applied to a spherical probe in 
medium with variable resistivity 
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0.00 0.01 £ AR )2 0.02 0.03 0.04 



where i* corresponds to R = 1.0 
i = i*, i* + 1 N-l , N 

The stability criterion for the cylindrical case is: 
M = 1 




where M is defined as in the spherical case as At/aR 2 . 

Computer program II (Appendix E) was used to generate the numerical 
solution to this problem. 
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IV. EXPERIMENT 



A. EQUIPMENT 

A radio frequency generator, Model RFG-2AV, manufactured by 
Radionics Incorporated of Burlington, Massachusetts, was used in the 
present investigation to generate radio frequency current. Output 
from the generator was a continuous, unmodulated sine wave at 500 
kilohertz with a maximum power of 18 watts. Although the generator 
included a built-in radio frequency voltmeter and radio frequency 
milliammeter , two independent RMS meters were utilized in data collec- 
tion to allow the use of more sensitive scales where required. Figure 2 
displays the generator, voltmeters, and one of the experimental probe 
test cells used in the experiments. Figure 3 is a schematic of the 
circuitry used. 

Both spherical and cylindrical radio frequency probe geometries 
were studied experimentally. An agar-water mixture (a gel-like 
substance) was used to simulate tissue. The mixture was contained in 
the test cells shown in Figure 4. In both cells, at least 25 probe 
radii separated the probe from the nearest solid boundary. Thus the 
cells could be considered as surfaces at infinity in relation to the 
probes. Plexiglass was used for the probe-agar cells due to its trans- 
parence and because it is a dielectric at the low power levels used 
during experimentation. Boundary conditions were changed from insulated 
surfaces to conductive surfaces by lining the inside surfaces of the 
cells with conductive foil. The current path was then from the probe, 
through the agar' medium, to the foil with electrical energy dissipated 
as heat in the agar. 
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Figure 2. Experimental Equipment 





Voltmeter V-j measured the voltage drop across the 
ten ohm resistor R-|q . Circuit current was deter- 
mined from the two quantities using Kirchoff's Law. 



Circuitry for radio frequency probe experiments 



Figure 3 
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Figure 4. Probe Test Cells 




Local temperature increases in the agar were measured by use of 
liquid crystals. Appendix B presents a complete discussion of liquid 
crystals. Briefly, liquid crystals change color at specific temperatures, 
the particular temperature depending on the crystal in use. The color 
change thus allows a continuous display of isotherm fronts for the map- 
ping of thermal regions. 

In the experiments, liquid crystals manufactured by National Cash 
Register and Hoffman-LaRoche were used. The Hoffman-LaRoche liquid 
crystal, ROCHROME, was provided in the form of a tape which could be 
applied directly to the object under study. The National Cash Register 
liquid crystals were obtained in liquid form. The crystals are actually 
encapsulated in spheres 5 to 30 microns in diameter, thus preventing 
contamination of the crystal. Due to this encapsulation procedure, it 
was possible to mix different liquid crystals together without changing 
the characteristics of any single crystal. Several distinct sets of 
isotherms were viewed since the liquid crystals mixed did not have tem- 
perature ranges that overlapped. 

In the experiments, three liquid crystals were mixed together and 
were applied to a transparent 0.003 inch mylar substrate. The liquid 
crystals were then covered with black paint and polyurathane varnish. 

The black background was required to absorb the light transmitted 
through the liquid crystals; the polyurathane varnish sealed the black 
paint. Sealing the system in this manner allowed repeated use of the 
liquid crystal system. Figure 5 displays the liquid crystal sheet 
encased in the holders used in the test cells. 

Prior to experimental trials, the crystals were calibrated using 
the Rosemount Constant Temperature Bath following the procedure of 
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Figure 5. Liquid Crystal Holders 






Petrovic [14]. The following liquid crystals were calibrated and 
subsequently used in experimental trials: 

Table 1 



Manufacturer 


Designation 


Red 


Green 


Blue 


National Cash Register 


R-27 


26. 7°C 


27.9°C 


28.9°C 


Hoffman-LaRoche 


Rochrome 


29. 9°C 


31 . 5°C 


34. 2°C 


National Cash Register 


R-37 


35. 7°C 


37.0°C 


38. 5°C 


National Cash Register 


R-53 


50.5°C 


52. 1 °C 


53.3°C 



Data were collected at specified times after application of power by 
photographing the isotherms displayed on the liquid crystal sheet. 

Distances from the probe center to an isotherm (color band) under consid- 
eration were easily measured from the photographs. 

B. EXPERIMENTAL PROCEDURE 

The agar solution was prepared by adding 3.0 grams of DIFCO 
"BACTO-AGAR" powder to 997 grams of boiling distilled water. When the 
powdered agar had thoroughly dissolved, a 100 ml. sample was taken to 
fill the resistivity test cell which was used to determine the variation 
of resistivity with temperature (see Appendix A). Prior to jelling, the 
remaining agar solution was poured into one of the probe test cells. The 
liquid crystal sheet was positioned and the probe located such that a probe 
diameter was in the plane of the liquid crystals. The entire apparatus was 
then cooled for approximately six hours to achieve a uniform temperature in 
the agar. At that time the initial agar temperature was recorded, the 
generator circuit was connected, and power was applied. Photographs of 
the liquid crystal sheet were taken through the agar at specified times, 
thus recording the color change in the liquid crystals. Distances were 
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measured directly from the pictures and v/ere scaled with respect to 
the probe radius to form a nondimensional distance. 

All quantities were nondimensional ized to allow general application 
of the data. The following nondimensional quantities were used (5): 



r . 



R 


_± 
»~ r 

0 


Distance 


T 

0 


at 

r z 
0 

T - T 

0 

T 

o 


Time 

Temperature 


G 


= fjY p o r o 2 


Power 


0 


V*7. kT 0 




'b 

P 


= 1 + me 


Resistivity 


where 


A 


= surface area 


of the probe (cm 2 ) 



for the sphere 4irr 2 

• o 

for the cylinder 2Trr L 

r. = distance from the probe center to the outer edge of the 
isotherm (color band) under consideration (cm) 
r Q = probe rad.ius (cm) 

'■ i 

a = thermal diffusivity (cm 2 /sec) 

t = time (sec) 

T = isotherm temperature (°C) 

T q = initial agar temperature (°C) 

I = radio frequency current (amps) 

p q = initial agar resistivity (ohm-cm) 
k = thermal conductivity (W/cm °K) 
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Using a constant temperature-power ratio (e/G^ ) , curves were 
plotted for distance (R) versus time (r) based on the solution of 
the finite difference equations for the probe geometry under investi- 
gation (Equations 8 or 11). The curves were then compared to the 
experimental data points of nondimensional isotherm location (R) at 
various times (x). An uncertainty analysis was conducted and is 
presented in Appendix D. 
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V. EXPERIMENTAL RESULTS 



A. SPHERICAL PROBE 

Figure 6 is a series of photographs of the growth of the radial 
temperature field around the spherical probe at nondimensional times 
of 1.5, 3.1, 6.2, and 12.5 (corresponding to times of 60, 120, 240, 
and 480 seconds for a spherical probe of radius 0.24 cm.). In 
Figure 6c, the isotherms used in data collection are labeled for 
reference. The temperature-color relationships used in data 
collection were: 

Liquid Crystal Temperature Color 



R-27 


27.9°C 


Green 


Rochrome 


29.9°C 


Red 


R-37 


37.0°C 


Green 



The R-53 appeared on the photograph at large times, but was not 
actually used in data reduction and comparison with theory. Values of 
nondimensional radius (R) and time (x) were plotted for comparison with 
the theoretical curves based on temperature-power ratios (e/G Q ). 

Figures 7 and 8 are representative of the comparisons. The dotted 
lines are the uncertainty limits, which are discussed in Appendix D. 
Experimental results were within 8 % of the theoretical values. Thus, the 
spherical theory, taking into account variable resistivity, was verified 
by the experimental procedure of using liquid crystals to determine 
thermal changes. 
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Figure 6a 

Tau =1.5 t = 60 sec. 



Figure 6b 

Tau =3.1 t = 120 sec. 




Figure 6. Spherical transient temperature field growth 
with probe diameter of 0.241 cm. 
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Figure 7 

Nondimens ional isotherm radii 
vs . 

nondimensional time for a spherical 
probe of diameter 0.482 cm. 
Data set I 




o 
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. 



Figure 8 

Nondimensional isotherm radii 
vs. 

nondimensional time for a spherical 
probe of diameter 0.482 cm 
using three liquid crystals. 
Data set J 
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B. CYLINDRICAL PROBE 

Figures 9, 10, and 11 display the experimentally determined 
transient temperature field at nondimensional times of 4.2, 6.25, 

9.3, and 12.5 (corresponding to times of 40, 60, 90, and 120 seconds 
for a cylindrical probe of radius 0.12 cm.). The length to diameter 
ratios used were 4, 10, and 20 in Figures 9, 10, and 11, respectively. 
Values of distance (R) were plotted as a function of time (t) for 
constant temperature-power ratios (e/G Q ), as seen in Figures 12, 13, 
and 14. The solid lines represent the theoretical curves. 
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Figure 9c 

Tau = 9.3 t = 90 sec. 



Figure 9d 

Tau = 12.5 t = 120 sec. 



Figure 9. Cylindrical transient temperature field 

growth, L/D = 4, probe diameter = 0.241 cm. 
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Figure 10a 

Tau = 4.2 t = 40 sec. 




Figure 10b 

Tau = 6.2 t = 60 sec. 





Figure 10c Figure lOd 

Tau = 9.3 t = 90 sec. Tau = 12.5 t = 120 sec. 



Figure 10. Cylindrical transient temperature field 

growth, L/D = 10, probe diameter = 0.241 cm. 



38 








Figure 11c 

Tau =9.3 t = 90 sec. 



Figure lid 

Tau = 12.5 t = 120 sec. 



Figure 11. Cylindrical transient temperature field 

growth, L/D = 20, probe diameter = 0.241 cm. 
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Figure 12 

Nondimensional isotherm radii 
vs. 

nondimensional time for a cylindrical 
probe with L/D = 4 and diameter = 0.241 cm. 
Data set M 




o 
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Figure 13 

Nondimensional isotherm radii 




o 
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Figure 14 

Nondimensional isotherm radii 
vs. 

nondimensional time for a cylindrical 
probe with L/D = 20 and diameter = 0.241 cm. 
Data set K 
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VI. DISCUSSION OF RESULTS 



A serious difficulty in modeling the radio frequency probe in 
tissue has previously been an inability to accurately map the thermal 
area surrounding the probe. Thermocouples and thermistors have either 
been placed inside the probe tip, where they actually measured the 
probe temperature, or placed in the medium, where they acted as thermal 
sinks. The use of liquid crystals, as demonstrated in this thesis, 
now allows an accurate, continuous display of specific temperatures 
depending on the crystal in use. Liquid crystals provide a flexible 
tool for temperature measurement regardless of probe geometry. Figure 15 
shows liquid crystals indicating isotherms which were caused by radio 
frequency current between two nonsymmetric probes. As can be seen, the 
liquid crystals react to provide a complete map of the thermal region. 

In the spherical geometry, the experimental values were in good 
agreement with the theory as modified for variable resistivity (8%). 

The cylindrical data and theory were not in agreement (30%). As a result 
of the poor agreement between theory and experiment in the cylindrical 
geometry, an examination of the flux lines emanating from the probe 
was conducted. Appendix C discusses the flux plots obtained. 

Figure 16 is a flux plot showing only the lines of constant current. 

As can be seen, the direction of the current lines is changed in the 
area above the probe. The bending of the current lines results in a 
higher local current density with greater heating than theoretically 
predicted. 
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Figure 15 

Liquid crystals map a 
nonsymmetric thermal region 
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Insulated Boundary 
(Top of Agar-water mixture) 




Figure 16. Constant current lines emanating from 
a cylindrical probe based on teledeltus plot of 
the isopotentials in plane perpendicular to probe 
axis 
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An approximation was derived to account for the increased current 
density. Due to symmetry of the problem, only half the probe-cell 
system was modeled. The following parameters were defined: 

Number of flow tubes (for half field) = b 
Current flow (for half field) = 1/2 
Half the probe surface area = ir^L 
Angle subtended by the ith flow tube = 

Surface area of the ith flow tube at the 
probe surface = <j>.r L 

Current flow per flow tube = I /2b = constant 
Current flux in the ith flow tube = I/2b<j>.r L 

1 O 

Recall that I was previously defined as 

I = , or for this situation, the local current density is 

t - I 
x i 2b*. r L 

1 O 



The local power parameter is G. *= 
for 1^, 



I. 2 p r 2 
1 0 0 

k T 



or after substitution 



G 1 " G o 

From Figure 16, * = 14°, b = 8 which yields G. = 2.58 G q . 

Applying this correction to the nondimensional temperature-power term 
( 9/G o ) shifts the theoretical curve so that experimental data more 
closely approximate the theory (Figure 17). 

Since the degree of current density is influenced by probe depth, 
the correction factor must be determined for each experiment. 

Examination of Figure 18 demonstrates the strong end effects on the 
finite cylindrical probe. The cylindrical theory (Equation 11) was 
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Figure 17 

Comparison of experimental data with 
original theoretical curve and 
corrected theoretical curve 

Cylindrical probe L/D = 20 

Data set K 
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Figure 18. Liquid crystals demonstrating the 
pronounced edge effects caused by a cylindrical 
probe (L'D = 20) 
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based on an infinite cylindrical probe with no end effects. The figures 
show the strong thermal fields at the probe ends developing before the 
effects at the probe center are observed. With time, the conduction of 
heat in the agar tended to reduce the thermal gradient between the probe 
ends and its midpoint, with the resultant growth as shown in Figure 11. 

A possible explanation for the early observation of thermal fields at 
the probe ends is that radio frequency current density is higher at 
sharp corners than at plane surfaces. The higher current density results 
in greater heating. 

By using a probe with a small surface area, resulting in high 
current density, it was possible to observe boiling in the agar. Since 
the liquid crystals used in the experiments were active in temperature 
ranges below the boiling temperature of the agar, no data were taken on 
boiling. However, Figure 19 is a photograph taken of the boiling 
phenomenon. It was of interest to observe formation of a bubble as the 
agar v/as heated to boiling. When rising, the bubble traveled until it 
encountered a cool region; the bubble then collapsed. This boiling is 
analogous to that occurring in actual surgery. The reason for controlling 
the probe tip temperature is readily apparent when observing the deteri- 
oration of the jelled agar as bubbles rise randomly towards the surface; 
the analogy in surgery being unpredictable physiological damage. 
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Figure 19. Boiling of the agar-water mixture 
caused by excessive probe tip temperature 
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VII. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

The use of liquid crystals for mapping regions surrounding heated 
surgical probes has been clearly demonstrated. When accuracy on the 
order of 1.0°C is sufficient, liquid crystals provide a continuous 
picture of thermal changes as displayed in Figures 20, 21, 22, and 23. 
Through the use of liquid crystals, the spherical theory developed by 
Gengler [5] has been verified. On the other hand, Gengler's cylindrical 
theory has been shown to require further work in order to account for the 
finite cylindrical probe with its strong end effects. 

It is most probable that an actual radio frequency surgical probe 
would be constructed with a length to diameter ratio in the range of 
two to five [1, 2, 15, 16]. Clearly, a theoretical treatment for an 
infinite cylinder would break down at these small length to diameter 
ratios. As a demonstration of this fact, an equivalent surface area 
for a spherical probe was determined based on the surface area of the 
cylindrical probe with a length to diameter ratio of four. Figure 24 
is a comparison of the cylindrical theory, the spherical theory, and 
the experimental data. Relationships between the cylindrical and the 
spherical parameters have been derived based on the fact that the cylin- 
drical surface area equals the spherical surface area ( 27rr C y] L= 4^>'5p^ ere ) • 

(6 ) i = I p o r cyl (G ) u = I p o r sphere 

o cyl k j 1 ■■ v o sphere 

0 0 
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Figure 20. Color photograph of liquid crystals in thermal 
field of cylindrical probe (L/D = 4). Current flux = 0.414 
amps/cm 2 




Figure 21. Color photograph of liquid crystals in thermal 
field of cylindrical probe (L/D = 10). Current flux = 0.414 
amps/cm? 
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Figure 22. Color photograph of liquid crystals in thermal 
field of cylindrical probe (L/D = 20). Current flux = 0.414 
amps/cm? 




Figure 23. Color photograph of liquid crystals in thermal 
field of spherical probe 
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Figure 24 

Comparison of cylindrical data (L/D = 4) and cylindrical 
theory with spherical theory and scaled cylindrical data 
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sphere 
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fore, 




L/2 from the equivalence of surface areas. 



There- 



^ G o^sphere ^o^cyl 
Similarly, 

R _ ^cyl _ T c,yl 

K sphere 2 T sphere 4 



As can be seen, the experimental data is more closely approximated by 
the spherical theory. 

The cylindrical probe with a length to diameter ratio of 20 most 
closely approximated the infinite probe of the cylindrical theory. 

Use of the correction factor to account for current path changes near 
the air-agar interface, tends to bring theory and experiment into 
better agreement. 



B. RECOMMENDATIONS 

Further work, both experimental and theoretical, is needed for the 
cylindrical geometry. Of particular interest are the cylindrical probes 
with length to diameter ratios of five or less since actual probes will 
be in this range. Experimental investigations into the effect of the 
air-agar interface with varying probe depths would demonstrate the 
analogous effects of lesion creation close to the skull; both the air- 
agar interface and the skull are essentially insulators. 

Liquid crystals placed perpendicular to the axis of the cylindrical 
probe and running from the air-agar surface to the bottom of the agar 
could effectively map the distorted thermal region due to the insulated 
air-agar surface. 



55 



Studies of various tissues should be conducted to fully catalogue 
tissue electrical properties. Such a catalogue would allow the surgeon 
to theoretically determine current strength and duration for a particu- 
lar lesion size prior to the actual surgery. 

Some refinement in the spherical experimentation would reduce 
uncertainty. The major uncertainties were in the thermophysical prop- 
erties of the agar. 

A method of retaining the liquid crystals in the plane of a diameter 
of the probe is required. This is particularly difficult when using the 
cylindrical probe. Placing the substrate in tension in its holder only 
succeeds temporarily. The substrate of mylar eventually stretched and 
sagged. Without having the liquid crystals in the diametral plane, true 
distances cannot be measured for scaling with respect to the probe radius. 

A better potentiometer is required for controlling the RFG-2AV gen- 
erator output. The control provided by the manufacturer had excessive 
play. This was partially overcome by careful selection of meter scales 
during experimentation. 

Use of a 35mm camera with a through the lense light meter system is 
suggested for taking data. Most data for this study was recorded using 
the industrial Polaroid film. When used under artificial light, the 
film is very slow (ASA 24); colors are not well defined which is criti- 
cal when using liquid crystals. The net result with the color Polaroid 
is a poor reproduction of the liquid crystal colors. Using 35mm black 
and white or tungsten corrected color film gives highly satisfactory 
results when used with the proper light meter settings. Cameras using 
Polaroid film are often moved during film insertion and removal. The 
slightest movement results in an unfocused picture at the short distances 
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normally used in photographing the probe area. With a 35mm camera, 
there is a minimum of movement during the shutter release and film 
winding processes. The actual time for inserting and removing the 
Polaroid film limits the data that can be taken; this is most true at 
small times. Using a 35mm camera, with the fast film advance time, 
allows more frequent data recording. 

A study should be conducted to determine the feasibility of mix- 
ing more than three liquid crystals without a loss of color brilliance. 
The required separation of temperature ranges should be determined 
since the closer the temperature ranges, the more crystals could 
theoretically be mixed. This would yield more data for each experi- 
mental trial. The temperature ranges, however, must be far enough 
apart so that each color band is displayed clearly. 
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APPENDIX A: RESISTIVITY 



BACKGROUND 

Initially, thermal and electrical properties of the 99.7 percent 
water/0.3 percent agar solution were unknown. It was decided to inves- 
tigate the electrical resistivity of the solution since thermal proper- 
ties could be correlated with the solution's water content [10]. A 
test cell, pictured in Figure 25, was designed and manufactured for 
use in the resistivity investigation. To prevent a standing wave in 
the resistivity test cell, the test cell length could not be an integral 
or a quarter wavelength of the 500 kilohertz signal used in the experi- 
ments. Using the relationship between frequency (f), wavelength (a), 
and the speed of light in air related to the speed of light in water 
(n) , 

3 x 10 8 

A = - meters 

n f 

the wavelength was determined and the test cell designed to prevent 
standing waves. 

A plexiglass tube,- 13 cm. in length with 2.5 cm. ID, was used as 
the test cell due to the insulating properties of plexiglass at the 
power levels used. Copper plugs were machined to fit tightly into the 
cylinder with the area of the plug ends known. A BLH two-inch, 
constantan-chromel needle thermocouple was introduced through the wall 
of the plexiglass cylinder midway between the copper ends. To prevent 
leakage of air into the container, a seal of neoprene rubber was provided 
between the test cell and the thermocouple shaft. During the experimental 
trials, the test cell was filled with hot agar solution. 
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Figure 25. The resistivity test cell with attached thermocouple 




By carefully capping the test cell, no air was trapped inside. 
The capped cylinder was placed in the generator circuit shown in 
Figure 26. 

Voltage was applied across the circuit for one minute when col- 
lecting data. One minute was chosen since no change in voltage drop 
occurred after this time. Determination of resistivity was accom- 
plished by recording voltages across the circuit and the ten ohm 
resistor along with the distance between the surfaces of the copper 
plugs. The voltage drop across the resistor yielded circuit current 
from Kirchoff's Law: 




The resistance of the cylinder was determined by dividing this current 
into the voltage drop across the cylinder. Resistivity is related to 
resistance by the area to length ratio: 



It was thus possible to calculate resistivity. 

At the same time the voltage readings were recorded, the needle 
thermocouple v/as read via the millivolt potentiometer. Data v/ere taken 
as the test cell cooled and again as the cell was heated in a controll- 
able oven (Figure 27). Early experimental trials showed a wide variation 
in agar resistivity (Figure 28). To reduce all plots to a common scale, 
two nondimensional quantities were introduced: 

T - T 

9 = j Temperature 

o 

p = Resistivity 
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MILLIVOLT POTENTIOMETER 




Voltmeter V-| measured voltage drop across the ten ohm resistor 
from which circuit current was determined. 

Voltmeter measured circuit voltage which was used in conjunction 
with circuit current to determine resistivity test cell resistance. 



Resistivity Test Cell Circuitry 
Figure 26 
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Figure 27. The resistivity test cell in the oven 
for the heating phase of resistivity data collection 
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Figure 28 

Variation in resistivity values 
as a function of temperature 
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where the subscript indicates the initial conditions of the experimental 
trial. In this way, the plots had common points of origin. 

A linear curve was fit to each set of data points. In dimensional 
form, 

p = p 0 + y (>.- V 

where y was the slope of the curve. Nondimensional izing yielded 

p = 1 -+ me 
yT 

where m = — - , the nondimensional slope. The values of m were 
p o 

generated by computer program III (Appendix E). Figure 29 is a sample 
of the nondimensional curves plotted by this program. 

As a result of this portion of the investigation, the great change 
in resistivity over the temperature range of interest (15-60°C) was 
demonstrated. It was determined that the programs written by Gengler [5] 
required modification to reflect the changing resistivity. This was 
accomplished by incorporating the slope of the linear curve from the 
resistivity data into computer programs I and II (Appendix E). 
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Figure 29 

Agar-water mixture nondimensional resistivity as a 
function of nondimensional heating and cooling temperatures 
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APPENDIX B: LIQUID CRYSTALS 



BACKGROUND 

Matter can exist in solid, liquid, or gaseous state depending upon 
the temperature and pressure of the surroundings. In the solid state, 
most compounds are composed of an ordered array of molecules in a 
crystal structure. When heated, the crystal begins to deform as the 
molecules absorb energy and leave the "ground state." The temperature 
at which the deformation begins is known as the melting point [17]; 
it is normally a well-defined characteristic of the compound. 

In 1888 Reinitzer noticed that not all compounds changed directly 
from solid state to an isotropic liquid [18]. The intermediate phase, 
a "mesophase," exhibited characteristics of both crystalline structure 
and isotropic liquid. Thus, the name of these substances was coined 
as liquid crystals [18]. Lehmann continued Reinitzer's work in liquid 
crystals. He classified the crystals into three categories based on 
molecular orientation; the categories are smectic, nematic, and 
cholesteric [17]. Only the cholesteric crystals will be discussed as 
it is the only type utilized in this report. Interested readers may 
find information on smectic and nematic crystals in the literature 
[16, 17, 18, 19, 20]. 

The cholesteric mesophase has the molecules' long axis parallel and 
coplanar. A model for this structure is a deck of cards, each card 
with one corner bent up such that a slight rotation of the entire 
stack results [18]. The average displacement of adjacent layers is 
fifteen minutes of arc per layer [21] with each layer on the order of 
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o 

1000 A thick [22], The rotation is due to the stretched out non- 
symmetric shape of the crystal molecule [23]. 

The helical rotation of the cholesteric liquid crystal gives rise 
to its optical properties. When transmitting polarized light, the 
helical layers cause the light to be rotated to the left with rotations 
of fifty degrees per millimeter. This property is known as being 
optically active. By way of comparison, quartz, one of the most 
optically active materials in nature, has a rotation of twenty degrees 
per millimeter for polarized light. Another property of the cholesteric 
liquid crystals is circular dichroism which is a function of temperature, 
material, and the angle of the incident light. The incident beam of 
unpolarized white light is split into two components. One component 
is rotated clockwise while the other component is rotated counter- 
clockwise. The result is that one component is reflected while the 
other is transmitted. To absorb the transmitted beam, a black background 
is required behind the liquid crystals. Circular dichroism is the char- 
acteristic that gives cholesteric liquid crystals their iridescent 
colors. This characteristic should not be confused with the birefringence 
or double defraction of white light by the cholesteric liquid crystals 
(Figure 30). Birefringence is often used as a test for the cholesteric 
mesophase [21]. In the material the incident beam is divided into two 
polarized components which are vibrating at right angles to each other 
and are travelling at different speeds. They are, therefore, refracted 
at different angles and emerge as two beams of perpendicularly polarized 
light. The rotary power of the birefringent cholesteric mesophase is 
as much as 1800 degrees per centimeter when light is transmitted along 
the optical axis normal to a thin sheet [23, 24]. The thickness of the 



67 




Birefringence: In the material, the incident beam of light is divided 
into two polarized components vibrating at right angles to each other 
and traveling at different speeds. They are, therefore, refracted at 
different angles and emerge as two beams of perpendicularly polarized 
light 




Circular dichroism: The incident beam of unpolarized light is split 

into two components, one rotating clockwise and the other rotating 
counterclockwise. The result is that one is reflected while the other 
is transmitted. 



Figure 30 



cholesteric film does not alter the predominant wavelength of the 
reflected light as it is rotated; however, the angle of incidence 
and the angle of observation both effect the wavelength of the reflected 
1 ight. 

The most striking changes in the cholesteric mesophase are those due 
to temperature changes. This mesophase is made temperature sensitive 
by the narrow band of optical reflection which is associated with a 
change in crystal structure. Cholesteric liquid crystals are also 
sensitive to organic compounds, shear, and electromagnetic radiation 
at high frequencies and power levels. Depending on composition, chol- 
esteric liquid crystals span a range from -20 to 250 degrees C [25]. 

Any specific compound may have its range narrowed to as little as one 
degree or increased to six degrees by varying the compound ingredients 
[26]. If heated far beyond its temperature range and allowed to cool, 
the transition colors occur at lower temperatures, similar to a magnet 
being heated above the Curie temperature [21]. If not thermally shocked 
as just described, the transition colors remain invariant with tempera- 
ture; the material may be repeatedly recycled. Not all cholesteric 
liquid crystals react with the same color changes; some display all colors 
of the spectrum starting with red to yellow to green to blue, while others 
may change directly from red to blue. 

Cholesteric liquid crystals will deteriorate in several days if not 
insulated from dust and organic materials. National Cash Register 
Company has been successful with encapsulating the cholesteric liquid 
crystals in spheres 5 to 30 microns in diameter [23]. The result of 
the microencapsulation process is a long life with little deterioration 
from contamination. The NCR process reduces the angle dependence of the 
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crystal as far as color is concerned although resolution is degraded 
and response time is increased approximately ten percent [22]. The 
microencapsulated crystals are used in industry since they are not 
readily contaminated; they were used in this experiment. 

Liquid crystals have many engineering applications. In temperature 
measurement, cholesteric liquid crystals are primarily used. Guidelines 
for use of liquid crystals in thermal mapping may be summarized as 
follows : 

a) The heat capacity of the object should be greater than that 
of the liquid crystals; 

b) The size of the object must be great enough to resolve the 
colors of the liquid crystals; 

c) The rate of temperature change must be slow enough for the 
liquid crystal to follow. (A maximum time of about 0.2 seconds.); 

d) The temperature range under study must be within the range of the 
1 iquid crystal ; 

e) The surface of the object must be capable of holding liquid 
crystals unless an adhesive backing is used; 

f) The surface must be black to display the liquid crystals; 

g) The temperature pattern must be correlated with another 
material property such as voids present in nondestructive 
testing or power dissipated per unit area. 
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APPLICATION AND CALIBRATION 

Reference 27 lists many techniques for application of liquid 
crystals; in this investigation, a commercial artist's air brush was 
utilized. Three crystals, designated R-27, R-37, and R-53, were mixed 
equally. After standing overnight, the cholesteric substances tended 
to separate from excess water. The water was carefully removed leaving 
a more concentrated liquid crystal which resulted in a more brilliant 
color display. The concentrated mixture was sprayed onto the transparent 
substrate in three thin coats, allowing adequate drying time for each 
coat. To prevent any absorption of moisture from the environment, dry- 
ing was conducted under a 40 watt lamp. After the three coats of liquid 
crystals were dried, two coats of black paint were applied. Testors 
Flat Black Spray Enamel was used since it contained no chemical which 
might have reacted with the crystals to cause a shift in the color- 
temperature relationship. The black paint was required to absorb the 
light transmitted through the liquid crystal , as only a small portion of 
the light was reflected. To provide a watertight system, a layer of 
polyurathane varnish was applied over the black paint and along each edge 
of the substrate. Thus, the crystals were sealed from all contaminants. 

A second type of liquid crystal , designated R0CHR0ME, was also used 
during the experiments. Rochrome, manufactured by Hoffman-LaRoche, 
consists of liquid crystals sandwiched between a layer of mylar and a 
black adhesive backing. Rochrome is cut to the desired shape and applied 
directly to the object under study. Again, the edges were sealed to pre- 
vent water contamination and separation of the Rochrome from the surface. 
The surfaces to which the Rochrome and the liquid crystals were applied 
in this investigation were mylar and glass. The mylar was used with the 
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probe apparatus while the glass was used during the calibration 
procedure. 

Calibration was conducted using the Rosemount Engineering Company 
Variable Temperature Bath (Model 913A) in conjunction with its 
commutating bridge and platinum resistance standard thermometer. 

National Cash Register liquid crystals (R-27, R-37, and R-53) 
were prepared by spraying them onto a glass slide, coating with 
black paint, and sealing with polyurathane varnish. Rochrome was 
prepared by sticking a small piece to the slide and sealing the edges 
with polyurathane. The slides were placed in the bath at approximately 
one-third of the bath's depth and close to the center of the bath 
diameter. Effects of the bath container edges were thus minimized. 

The bath temperature was raised slowly as the slides were observed. 

To minimize variations between observers, the same observer read the 
slide colors while an assistant operated the equipment during all 
calibration procedures. 

As the temperature was raised, the liquid crystals reacted, each 
at its own temperature level. When the slide had turned one color, 
i.e. red to green, the commutating bridge reading was recorded; then 
the temperature of the bath was increased to determine the next color 
change. A conversion table was used to convert the bridge readings to 
temperature. The uncertainty in this temperature-color calibration 
procedure was calculated to be ±0.5°C which agrees with Reference 27. 
Table B-l lists the crystals used and their calibration temperature- 
color relationships. 
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Table B-l 



Crystal 


Red 


R-27 


26.7°C 


Rochrome 


29.9°C 


R-37 


35.7°C 


R-53 


50.5°C 



Green 


Blue 


27.9°C 


28. 9°C 


31 . 5°C 


34. 2°C 


37.0°C 


38. 5°C 


52.1°C 


53.3°C 
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APPENDIX C: FLUX PLOTS 



Gengler's theoretical analysis of the cylindrical probe modeled the 
probe as an infinite cylinder with negligible end effects; the current 
flow should have been a one-diniensional problem. With the experimental 
cylindrical probe, end effects were pronounced, prompting an investiga- 
tion of current paths in the probe vicinity. The area of interest was 
within ten probe diameters. The problem was approached in two phases. 
First, the horizontal plane was examined using TRUMP [28, 29], a computer 
code, while in the second phase, teledeltus paper was used to depict the 
isopotentials in the vertical plane. 

TRUMP solved the general partial differential equation v 2 e = 0. 

It was postulated that the radio frequency potential steady state 
condition was established instantaneously with power application. 

A nodal network was constructed’ for use with TRUMP (Figure 31). Due 
to symmetry of the probe-cell system, only one quadrant of the system 
was modeled. Three cylindrical probes, with length to diameter ratios 
of 4, 10, and 20 were used with two sets of boundary conditions. The 
first boundary conditions approximated a small finite cylinder inside 
an infinite cylinder by the use of two conductive sides of the nodal 
boundaries for the outer cylinder wall and two insulated sides of the 
nodal boundaries for the endcaps of the outer cylinder. Figures 32, 

33, and 34 are plots of the isopotentials and lines of constant current 
for each probe length to diameter ratio. Inspection of these figures 
in the area within ten probe diameters, shows the current lines orthogonal 
to the probe surface. It was, therefore, expected that the current flow 
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Probe (L/D = 20) 
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TRUMP nodal network 
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and the resultant thermal problem would be one dimensional. For 
comparison, a second set of boundary conditions surrounded the probe 
with conductive surfaces. Figures 35, 36, and 37 display the iso- 
potentials and constant current lines for this situation. 

The second portion of current path investigation was undertaken 
when experimental data did not agree with the theory. As stated earlier, 
in this second phase, isopotentials in the vertical plane were examined 
using the same sets of boundary conditions as in the first phase. 

The probe-cell system was scaled on teledeltus paper. Boundary 
conditions were simulated using conductive silver paint and metallic 
tape. Figures 38 and 39 are plots of the isopotentials and constant 
current lines under the two sets of boundary conditions. It is seen 
that the insulated agar surface caused current to bend towards the 
conductive boundary. Since heating is a function of current per unit 
area, the higher current density resulted in higher heating effects 
in the experimental model. The "higher heating effects were the major 
cause for the difference between theory and experiment with the 
cylindrical probe. 
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Figure 35 
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<f>‘ = 0 (Free surface of agar-water mixture) 




Figure 38 

Teledeltus plot of 
isopotentials and constant 
current lines in a plane 
perpendicular to the cyl- 
indrical probe axis. 

Free surface of agar- 
water mixture acting as an 
insulator. 



(Free surface of agar-water mixture) 



4 > = 0 



Figure 39 

Teledeltus plot of 
isopotentials and con- 
stant current lines in 
a plane perpendicular 
to the cylindrical 
probe axis. 

<f>'=C 

Free surface of agar- 
water mixture acting 
as conductive boundary. 




<f> = 0 
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APPENDIX D: SAMPLE UNCERTAINTY ANALYSIS 



Reference 30 categorizes error in experimental investigation as 
fixed error, variable error, and mistakes. Fixed error is defined 
as a constant error such as equipment consistently reading low. 
Variable error is caused by variations in sensors, an example of 
which is having two different observers taking readings in separate 
experimental trials. Finally, mistakes may occur during experiments 
by reading the wrong scale, etc. If it is assumed that the error, 
i.e. uncertainty, of each primary variable has the same probability 
of occurrence, it is then possible to calculate the overall uncer- 
tainty in each value. As an example, define P, the result of some 
experiment as, 



uncertainty in each primary independent variable, then the relative 
uncertainty in P is defined as 



In the experimental trials on the spherical probe, the quantities 
calculated were 



P = P(x-j, x 2 , x 3 ,...x n ) 



If Wp is the uncertainty in P, and u-j , . ..w n are the 
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Deriving the expressions for the relative uncertainty in each 
quantity 




For all independent variables except resistivity p Q it was possible 
to estimate the uncertainty either directly from observing equipment 
(voltage, current) or from available tables (a, k). For resistivity 
it was necessary to calculate uncertainty using 




Following the example of Equation D-l , 
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Table D-l lists values of uncertainty for each independent variable 
and the value of the variable used in the calculations. The following 
values were calculated for uncertainty: 



03 




R 



03 

T 

T 



0.036 (3.6%) 



0.1217 (12.2%) 
0.019 (1.9%) 
0.029 (2.9%) 
0.102 ( 10 . 2 %) 



To obtain 
was used 



uncertainty in the quantity 6 /G q the following relationship 




Using the values above. 



°Vg 

-eTe 2 - = 0.123 (12.3*) 

0 
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As may be seen from the last column of Table D-l, the largest 
uncertainties arise from the thermophysical properties, i.e. thermal 
conductivity (k), thermal diffusivity (a), and resistivity ( p o). 

Table D-l 

Y / \ 

Quanti ty i Uncertainty i Percent Uncertainty 



V cct 


= 


80.0 volts 


0.0025 volts 


0 


V R 


= 


0.118 volts 


0.0025 volts 


0.02 


D 


= 


2.540 cm 


0.025 cm 


0 


L 


= 


11.95 cm 


0.025 cm 


0 


R 10 


= 


10.0 ohms 


0.3 ohms 


0.03 


V 


= 


0.60 volts 


0.0025 volts 


0 


r 

0 


= 


0.241 cm 


0.001 cm 


0 


T 

0 


= 


23.0 °C 


0.25 °C 


0.01 


T 


= 


37.0 °C 


0.10 °C 


0 


r 


= 


0.380 cm 


0.01 cm 


0.02 


° P 

r p 


= 


0.90 cm 


0.01 cm 


0.01 


t 


= 


30.0 sec 


0.025 sec 


0 


k 


= 


0.628 W/m°K 


0.0628 W/m°K 


10.0 


a 


= 


0.001512 m 2 /sec 


0.0001512 m 2 /sec 


10.0 


P o 


= 


3250.0 ohm-cm 


117.0 ohm-cm 


3.60 
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APPENDIX E 



C COMPUTER PROGRAM I 

C PROGRAM FOR THE INSULATED SPHERE. 

C 

C 

C MODIFIED TO REFLECT CHANGING RESISTIVITY 

C 

C 

c 

C TET REPRES'ENTS THE NONDIMENS IONAL TEMPERATURE. 

C BETA IS THE BLOOD FLOW TERM. 

C G IS THE NONDIMENSIONAL POWER TERM. 

C DR IS DELTA RADIUS. 

C EM REPRESENTS THE M IN THE STABILITY CRITERIA. 

C DELTAU REPRESENTS NONDIMENSIONAL TIME CHANGE. 

DIMENSION TFT( 500.2 ) ,TT(500.2) t RH0(500) 

I T I M E = 0 

RE AD (5 , 100 ) BET A, DR, GO. AM, L I , IH 
100 FORMAT! 3F10. 4, F10. 6,21 10) 

WRITE(6,150)BETA,DR,GG,AM,LI,IH 
150 FORMAT! • 1 • ,6X BETA = ‘ , F4 . 2, / , 7X» 1 DR = • , F6 . 4, / , 7 X, • G 

AO' = « ,F6. 1 ,/,7X, • SMALL M = • , F 1 0 . 6 , / , 7X . * L 1 = *,I3,/,7X 
B , 1 HI = * , 1.4* / / ) 

TETtl ,1 ) =0 . 0 
WRITE! 6,2000) 

2000 FORMAT !• 1» ,T4, ‘TEMPERATURE DISTRIBUTION AT SPECIFIED T 
A I ME AND RADIAL L OC AT IONS • , / / , 9 X , * T IME ( T AU ) • , IX , • LOCAT 
B I ON ( R 1‘, 3X. • TEMPERATURE ( THETA ) /POWER (GO)*) 
S=AM*DR**2*G0 
EM=1. 0/(2. 0+BETA ? :-DR*F2-S) 

DELTAU* FM--‘ D 5 ' ■ 2 
DO 10 K = 1..IH 
10 TET! K , 1 } =0 . 0 

50 DO 200 I=L I. IH 

TET! IH+-1.1)=TET{ IH-1.1) 

TFT! LI-I , L)=TETt LI + 1 ,1) 

X=FLGAT ( I ) 

200 TETtl .2)*EM* ( ( l • — 1 . /X ) *TET ( 1-1 , 1 ) + ( 1 . +1 . / X )4 ! T ET ( I + 1* 1 ) 

fi ) +TET { I • 1 )v! l.-2.*EM+< EM&DR*&2* AM*GO) / ( X*DR) *«4) + ( GO-"- E 
CM* DR* *2 >/ ! X' DR )*-*4 * 

I TIMF=T TIMF+I 
T I M E = FLOAT ( I TIME ) 

T AU=1 I M L ; r DEI TAU 
IF! ITIME.GT.1000) GO TO 90 
I F (IT IMF. GT. 100) GO TO 70 
■I TEST = I II ME/20 F 20 
GO TO 80 

70 I T EST = I T I M E/ 50^5 0 
GO TO 80 

90 ITFST=I TIME / 1 00v< 100 

80 J = 0 

DO 60 J=LI , IH 
Y=FLOAT ( J ) 

R=Y*DR 

IF ! < I TI ME- 1 TEST ) . NE.O ) GO TG 60 

TT( J , 2 ) =TE T( J , 2 ) /GO 

WRITE ( 6 ♦ 220 ) TAU , R, TT < J , 2) 

220 FORMAT!' «. , F 1 5 . 4 , 7X , F6 . 3 , 5 X , El 5 . 7.) 

60 TET! J , I ) = TET( J, 2) 

I F( IT IMc.GT.4000 )STOP 
GO TO 50 
END 
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C COMPUTER PROGRAM I I 

C PROGRAM FOR THE INSULATED CYLINDER 

C 

C MODIFIED FOR VARIABLE RESISTIVITY 

C 

C TFT REPRESENTS THE NONDIMENS IGNAL TEMPERATURE. 

C BETA IS THE BLOOD FLOW TERM. 

C G IS THE NONDIMENS IONAL POWER TERM. 

C DR IS DELTA RADIUS. 

C EM REPRESENTS THE M IN THE STABILITY CRITFRIA. 

C DELTAU REPRESENTS NOND I MENS IONA L TIME CHANGE. 

DIMENSION TFT I 500. 2) ,TT(500,2) 

IT IME = 0 

READ! 5 . 100) BETA, DR, GO. AM, LI . 1H 
100 FORMAT ( 3F10.4.F 10. 6. 21 10) 

WRITE (6. 150 ) BETA, DR. GO, AM, L I . IH 
150 FORMAT ! ' 1 6X » ' BETA = ' , F 4 . 2 , / , 7X , ' DR = ' » F6. 4, / , 7X, ' G 

AO = ',F6.3./,7X. ' SMALL M = • ,F 1 0. 6 , / , 7X , • L I = ',13,/,' 
BHI = • , 14 , // ) 

WRITE! 6,2000 J 

2000 FORMAT! • IE, 14. 'TEMPERATURE DISTRIBUTION AT SPECIFIED 
A I ME AND RADIAL LOCAT IONS • . // .9X, • T IME ( T AU ) ' , IX , • LOC 
B I ON ! R 3X. • TEMPERATURE ! THETA )/ POWER (GO)') 

EM= I .0 / ( 2 . 0-8 ET A&D R4#2-AM £ GO ) 

DE LTAU = EM'? DR>*2 
DO 10 K=1 . IH 
10 T FT !K. I ) = 0 .0 
50 DO 200 I I . IH 

TFT! IH+I. I) = TET( I H- 1 , 1 ) 

TET !LI—1,L)=TET ! L I +1 » 1 ) 

X= I 

200 TET! I, 2) = EM-: ! ( I . - 1 . /( 2 . X ) ) * TE T ( I -1 , 1 ! + ( 1 . + 1 . / ( 2 . * X ) ) •:= 
ATET ( I +1 , 1 ) )+< I .-BETA-, E M - D R ^ * 2 - 2 . T E M +( A M ■> GC ; ; EM > / X-YW2 ) vT 
BET! 1 . 1 1 +( GO" E MJ / X'££2 
I T IME = I T I M E+ I 
T AU = I T I MF'T DELI AU 
I F! IT IME. GT. 1000) GO TO 90 
I F ( IT I M E. GT . 100 i GO TO 70 
I TEST = I T I M F / 2 0 2 0 
GO TO 80 

70 I T EST =ITI ME/50 *50 

GO TO SO 

90 ITEST=ITIME/100*100 

80 J =0 

DO 60 J=LI , IH 
R = J 4 D R 

IF! ! IT IME-ITEST ) .NE.O ) GO TO 60 
TT< J.2) =TFT! J,2) /GO 
WRITE(6,220)TAU.R.TT(J. 2) 

220 FORMAT!' ' , FI 5 . 4 , 7X , F6 . 3 . 5X , El 5 . 7 I 
60 TFT! J.1)=TET! J,2) 

IF! IT IME. GT. 4000 JSTOP 
GO TO 50 
END 



89 






onnnnnnooooo 



PROGRAM NUMBER III 

PROGRAM TO FIT AND PLOT CURVES FOR RESISTIVITY AND 
TEMPERATURE DATA. SUBROUTINE LS0PL2 IS CALLED TO FIT 
THE CURVES TO DATA. SUBROUTINE DRAW IS USED TO PLOT 
THESE CURVES AND THE DATA POINTS. ALL ARGUMENTS FOR 
LSOPL2 MUST BEGIN WITH D EXCEPT M AND KM. 

FOR ALL DRAW THE FOLLOWING NOTATION IS USED 
M=NUMP T S . THE NUMBER OF POINTS 
T EMP= NOND IMFNS I ONA L TEMPERATURE (X) 

RES = RES 1ST I V I TY (Y) 

DX REPRESENTS NONO I MENSI ONAL TEMPERATURE— -THETA . 

DF2 IS RESISTIVITY 
IMPLICIT REAL*8(D) 

REAL* 8 TITLE! 10) /1 0'*' •/ 

REAL-1‘8 ITITLE 
REALS8 LAB, EL/' •/ 

REAL LABEL 2/4HHEAT/ 

REAL LA3FL1/4HCGGL/ 

DIMENSION D X ( 100) , DF2(100) , DWI ( 100) , DY ( 1 00 ) .DELY(IOO), 
ADB(IOO) i DS B( 100) «THET( 100). TEMPI 100), RES! 100) , B ( 1 0 0 ) »Y 
B I 100) . I TITLE 1 12) , DTHET ( 103 ), DRES ( 100 ), THET A ! 100 ) 

READ! 5. 9)1 I TITLE I I ) ,1 = 1,12) 

9 FORMAT (A3) 

no ioo j=i . 7 

READ! 5. 1 >M .KM. XCAL , YCAL 

1 FORMAT I2I10.2F13 .L) 

WRITE! 6,2) M. KM, YCAL, XCAL 

2 FORMAT ( 30HI THE NUMBER OF DATA POINTS IS .I10./.30H0 T 

AHE MAXIMUM DEGREE OF FIT IS . I10./.21H0 THE Y-AXIS SCA 
BLE IS .F10.1.21H0 IN DEGREES PER INCH ,/,21H0 THE X-AX 
CIS SCALE IS . E L 0 . 1 , 20HO IN OHM-CM PER INCH ) 

RF AD ( 5 « 3 )( DX ( I) , DF2 ( I ) , I = 1 , M ) 

3 FORMAT! 2F10. I) 

R E AD! 5 . II )T0 

11 FORMA T ( F 5* l ) 

WRITE I 6. I2)TG 

12 FORMAT ( 34H0 THE INITIAL TEMPERATURE (TO) WAS, F5. I , 9H0I 
D TG=TQ 

DO 7 1=1. M 

ORES! I ) =QF2 ( I )/DF2 ( I ) 

DTHET! I ) = ( OX (I )-DTO)-/DTO 

7 DWI (I)=i. 0000000 
WRITE (6.4) 

4 FORMAT { 14X.45H0DIM5NSI0NAL VALUES NGNDIMENSIONAL VA 

BLUES ./.9X.51H TEMPERATURE RESISTIVITY THETA 

C RHD-PRIMF. , / . 9X . 24H (DEGREES C) (OHM-CM) ,//) 

WRITE! 6.6) (DX! I ) .DF2! I ) .DTHET! I ) .ORES! I ) , I = 1 , M ) 

6 FORMAT! 12 X , F6 . 1 . 7X , F6 . 1 , 8X . F8 . o , 5X.F8.6) 

CALL L SOP L2(M, KM, DTHET, DRES, DWI , DY , DEL Y , D8 , CS B , T I T L E ) 
DO 15 1=1, KM 
15 B(I)=D3(I) • 

DO 8 I = 1 * M 
Y ( I ) = DY ! I I 
THETA U )= DTHET ( I ) 

8 RES! I ) =ORES( I ) 

IF!J.EQ.2) GO TO 101 

99 CALL DRAW CM. THETA, RES, 1,5, LABEL, ITITLE, XCAL,YCAL, 0,0,0 
A, 0. 8, 1 1 . 1..LAST) 

CALL DR AW I M. THET A, Y, 2,0, LAB EL 2, ITITLE, XCAL, YCAL,0, 0,0, 
AO, 8,11, 1 . LAST) 

WRITE(6« 13) LAST 

13 FORMAT ( 8 HO LAST = ,13) 

GO TO ICO 

101 CALL DRAW! M.THETA, RES, 2, 4, LABEL, ITITLE, XCAL, YCAL, 0,0,0 
A, 0,8. 11 ,1,-LAST) 

CALL DRAW! M.THETA. Y,3, 0, LAB ELI ,1 TI TLE , XC AL , YCAL , 0 ,0,0, 
AO, 8. 11, I. LAST ) 

WRITE (6.1 04) LAST 
104 FORMAT! 8H0 LAST = ,13) 

100 CONTINUE 
STOP 

END 
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APPENDIX F 



TEMPERATURE AND 
DATA 



RESISTIVITY 
SET A 



VALUES 



HEATING CONDITIONS 



THE NUMBER OF DATA POINTS IS 5 

THE MAXIMUM DEGREE OF FIT IS 1 

THE Y-AXIS SCALE IS 0.10 IN DEGREES PER INCH 

THE X-AXIS SCALE IS 1.00 IN OHM-CM PER INCH 

THE INITIAL TEMPERATURE (TO) WAS 21.00IN DEG C 



DIMENSIONAL VALUES 
TEMPERATURE RESISTIVITY 
(DEGREES C) (OHM-CM) 



27.5 
31.3 

37.5 
44.0 
51'. 4 



2753.6 
2508.5 
2508.5 
2014. 9 
1800.1 



NONDIMENS IONAL VALUES 
THETA RHO- PRIME 



0.309524 
0.490476 
0 .785714 
1.095238 
1.447619 



1 .000000 
0.910989 
U. 910989 
0.731733 
0.653726 



COEFFICIENTS OF THE POWER SERIES EXPANSION 
Y(X)=B(1)+6(2)*X+B(3 )*X**2. . . 

R ( 11= 1.09146D 00 B ( 2 ) = -3. 02738D-01 

ESTIMATES OF ERROR FOR THE COEFFICIENTS 



E RRB ( 11 = 9.2330-02 E R R 6 ( 
SUM SQ DEV = 1 . 0233—03 
CHI SQ= 5.993D-03 

I X( I) F2( I ) 

1 0.310 1% OOOOOOOD 00 

2 0.490 9 .1098925D— 01 

3 0.786 9. 10989250-01 

4 1.095 7.317 3300D-0 1 

5 1.448 6>. 5372603 D-01 



2 ) = 1 . 008D-0 1 

E- RAT IC 
DEG OF FREEDOM 

Y ( I ) 

9.977581 70— 01 
9. 4297694 D-01 
3.5359704 D— 01 
7. 59 8 923 ID- 0 1 
6 .53 213 0 70—01 



; 9.028D 00 

3. 

DELY ( I ) 

•2. 241 8303D-03 
3.19 8 769 ID- 02 
-5 . 73 922 OS 0-02 
2.81 59305D-02 
•5 . 1295695D-04 
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TEMPERATURE AND RESISTIVITY VALUES 
DATA SET A 

COOLING CONDITIONS 



THE” NUMBER OF DATA POINTS IS 9 

THE MAXIMUM DEGREE OF FIT IS 1 

THE Y-AXIS SCALE IS 0.10 IN DEGREES PER INCH 

THE X-AXTS SCALE IS 1.00 IN OHM-CM PER INCH 

THE INITIAL TEMPERATURE (TO) WAS 21.00IN DEG C 



DIMENSI 
TEMPERATURE 
(DEGREES C) 

25.2 
29:. 9 

31.2 

32.7 
33*. 6 

34.6 
35.9 

40.8 

44.7 



ONAL VALUES 
RESISTIVITY 
(OHM- CM) 

2870.5 

2687.9 

2645. 9 
2546.3 

2508.5 
2453.8 

2418.7 
2242 .2 

2102.7 



NON DIMENSIONAL VALUES 
THETA RHC-PRI ME 



0.200000 
0.423810 
0.485714 
0.557143 
0 .600000 
0 . 6476 1 9 
0.709524 
0 .942857 
1.128571 



1.000000 
0.936387 
0.921756 
0.887058 
0.873390 
0.854834 
0. 842606 
0.781113 
0.732520 



COEFFICIENTS OF THE POWER SERIES EXPANSION 
Y ( X ) = 3 ( 1 ) +B ( 2 ) 1 X + B ( 3 ) * X^*2 . . . 

R( 1)= 1.05463D 00 B( 2) =-2. 9 1 728D-0 1 

ESTIMATES OF ERROR FOR THE COEFFICIENTS 



F RRB ( 1 ) = 1 . 928D- 02 ERRB( 2) = 2.318D-02 

SUM SO DEV = 4.1 700-05 F-RATIO = 1.071D 02 

CHI SQ= 4 . 3520— 04 DEG OF FREEDOM = 7. 



1 
1 

2 

3 

4 

5 

6 

7 

8 
9 



X(I) F? ( 1 ) 

0.200 1 . OOOOOOOD 00 

0.424 9 . 36 38 l 39D- 0 1 

0.486 9 .21 75579D-01 

0.557 8.87C5800D-01 

0.600 8. 73889570— 01 

0.648 8. 5483365D— Oi 

0.710 8^ 42 605 8 2D— 01 

0. 943 7 .-81116270-01 

1. 129 7. 3252047D-0L 



Y ( I ) 

9. 96279840-01 
9. 309883 6D-01 
9.129290 ID— 01 
8. 92 09131 .0-01 
6. 7958668 D— 01 
L .6569688 C-01 
8 . 47637 5 3C— 01 
7 . 79 567 69D-01 
7.25 38965 D— 0 1 



DEL Y< I ) 

-3 . 72 01594 D— 03 
-5. 3990293D-03 
-8. 82677790-03 
5 . 0333033D-03 
5. 69911680-03 
1 .08632260-02 
5.03 171 5QD-Q3 
-1. 55058110-03 
-7 . 1308133D-03 
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TEMPFRATJRE AND RESISTIVITY VALUES 
DATA SET B 



HEATING CONDITIONS 



THE NUMBER OF DATA POINTS IS 5 

THE MAXIMUM DEGREE OF FIT IS 1 

THE Y-AXIS SCALF IS 0.10 IN DEGREES PER INCH 

THE X-AXIS SCALE IS 1.00 IN OHM-CM PER INCH 

THE INITIAL TEMPERATURE (TO) WAS 23.00IN DEG C 



DU ME NS l 
TEMPERATURE 
( OEGRF'ES C) 

26.2 
32-. 2 
AO . 4 

46.0 

52.0 



ONAL VALUES 
RESISTIVITY 
(OHM-CM ) 

3723.4 
3422.1 

2997.6 

2753.6 

2500.4 



NONOI MENS IONAL VALUES 
THETA RHO-PRIME 



0. 139130 
0.400000 
0.756522 
1 .000000 
1 .260870 



1.000000 
0.919079 
0.805071 
0. 739539 
0.673685 



COEFFICIENTS OF THE POWER SERIES EXPANSION 
Y( X) =8(1 )+6(2) *X+B(3 J ‘ X‘ *2. . . 



B ( 1)= 1.03617D 00 B ( 2 ) =-2. 93393D-0 1 

ESTIMATES GF ERROR FOR THE COEFFICIENTS 

F RRB ( 3 ) = 1 .53 3'D— 02 ERRB( 2)=1.875D-02 
SUM SO JcV = 3.4240-05 F-RATIO = 2.449D 02 



CHISQ= 2 


.2180-04 


DEG OF FREEDOM 


= 3. 


I 


X( I I 


F2( I ) 


Y ( I ) 


DELY ( I ) 


1 


0.139 


1 . 000000 OD 00 


9. 9 5 34 6 5 4D— 01 


-4. 65346390-03 


2 


0.400 


9 s. 1907934 D- 01 


9.1880931 D-01 


-2.70025330-04 


3 


0.7 57 


8. 0507063D— 01 


8. 14208440-01 


9.13730180-03 


4 


1.000 


7 . 395391 3D— 01 


7. 4277369 D-01 


3. 23456160-03 


5 


1. 261 


6* 7368534 D-01 


6.66 23647 D-0 1 


-7.44887420-03 



93 



TEMPFRATURE AND RESISTIVITY VALUES 
DATA SET B 



COOLING CONDITIONS 

THE NUMBER OF DATA POINTS IS 19 

THE MAXIMUM DEGREE OF FIT IS 1 

THE Y-AXIS SCALE IS 0.10 IN DEGREES PER INCH 

THE X-AXIS SCALE IS 1.00 IN OHM-CM PER INCH 

THE INITIAL TEMPERATURE (TO) WAS 23.00IN DEG C 



DIMENSIONAL VALUES 


nondimens ional 


VALUES 


PERATURE 
GREES C) 


REST STI VI TY 
(OHM-CM ) 


THETA 


RHO-PRIME 


24.1 


3872.4 


0.047826 


1.000000 


24 .7 


3850.4 


0.073913 


0. 994319 


25.7 


3764.8 


0.117391 


0.972214 


26.5 


3743.9 


0.152174 


0.966816 


2B. 3 


3662.9 


0.230435 


0.945899 


20.4 


3585.3 


0.278261 


0.925860 


30.0 


3547.7 


0.304348 


0. 916150 


31.0 


3492.8 


0 .347326 


0.901973 


32.5 


3394.7 


0.413043 


0.8 76o40 


35.7 


3226.3 


0.552174 


0. 833153 


3S-.0 


3136.6 


0 .652174 


0.809989 


40.3 


3079.5 


0.752174 


0.795243 


42 .0 


2971.3 


0.826087 


0. 767302 


46.2 


2822 .5 


1 .00e&96 


0.728876 


48.0 


2753.6 


1.086957 


0.711064 


51.4 


2546.2 


1.234783 


0 . 6 5 76 2.0 


53.5 


2527.2 


i .326087 


0.652619 


56.3 


2453.0 


1 .447826 


0. 633664 


5 9.6 


2384.6 


1 .591304 


0.6157S4 



COEFF ICT ENT S OF THE POWER SERIES EXPANSION 
Y( X) =B ( l)+8( 2) :'-X+B(3) *X**2. . . 

R( 1 )= 9 . 96 85 3D— 0 L 8< 2 ) =-2 . 59 982D-0 1 



F R RB ( 


ESTIMATES OF ERROR 
l ) =2 . 4G4D-02 E RRB ( 


FOR THE COEFFICIENTS 
21=2.9470-02 


SUM 


SO DEV 


= l. 8400-04 


F-RATIO 


= 7.784D 01 


CHISO= 4 


.7870-03 


DEG OF FREEDOM 


= 17. 


I 


X ( I ) 


F2( I 1 


Y ( I ) 


DELY(I ) 


1 


0.048 


1 . OOOOOOOD 00 


9.3441862D-01 


-1. 55813760-02 


2 


0.074 


9. 94318770-01 


9.77 63649 D-01 


-1 .66822 34D-02 


3 


0-117 


9.722 136 1 D-01 


9. 6 63 32 9 2 C— 01 


-5. 88 06934D- 03 


4 


0.152 


9 .66816440-01 


9. 57290070-01 


-9 . 5263750D- 03 


5 


0.230 


9 . 45 6991 3D— 0 1 


9. 36 943 6 5 D-01 


-8 . 95 5 52 99 D— 03 


6 


0.278 


9. 2585993D— 01 


9. 2450973D— 01 


-1. 3501984D-03 


7 


0.304 


9. 1615019 D— 01 


9.1772759 D— 0 1 


1 .57740390-03 


8 


0.343 


9. 01 9729 4D- 01 


9. 0642403D-01 


4.45109410-03 


9 


0.413 


8. -766393 1 D— 0 1 


8 . 6 9 4686 8 D— 01 


1 . 28 2 86 75D- 02 


1 0 


0.552 


8.33152 57D— 01 


8.53 29728 D-01 


2.01 447 12D-02 


1 1 


0. 652 


8.:0958864D-01 


8. 2729908D— 01 


1. 73 10444D-02 


12 


0.75 2 


7 .'952 43 2 6 D— 01 


8.0 1300860^01 


6.057 6240D- 03 


1 3 


0.826 


7 . 673019 3D — 01 


7. 82 0848 2 D-01 


1 . 47828 93 D— 0 2 


14 


1.009 


7u 288761 5 D-01 


7. 346098 6D- 01 


5. 7337062D-03 


1 5 


1.087 


7. 11 083570-01 


7. 14263440-01 


3 . 1798741D-03 


16 


1.235 


6. 5752505D-01 


6. 758313 2D-01 


1 . 6306273D-02 


17 


1.326 


6. 5261853D-01 


6. 520936 4D— 01 


-5. 246942 3D-04 


1 8 


1 . 44 3 


6. 3366388D-01 


6 . 204438 6D-01 


-1 .32200200-02 


1 9 


1.591 


6. 157 9382D-0 1 


5.331421 0D-01 


-3. 2651727D-02 
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TEMPERATURE AMD RESISTIVITY VALUES 
DATA SET C 



HEATING CONDITIONS 



THE NUMBER OF DATA POINTS IS 5 

THE MAXIMUM DEGREE OF FIT IS 1 

THE Y-AXIS SCALE IS O.IO IN DEGREES PER INCH 

THE X-AXIS SCALF IS I. 00 IN OHM-CM PER INCH 

THE INITIAL TEMPERATURE (TO) WAS 22.50IN DEG C 



DIMENS I 
TEMPERATURE 
(DEGREES C) 

24.8 
29.5 
35 .4 
41.3 

47.8 



GNAL VALUES 
RESISTIVITY 
(CHM-CM) 

3428.4 

3110.6 

2799.1 

2544.3 

2299.9 



NONDIMcNSIONAL VALUES 
THETA RHO-PRIME 



0. 102222 
0.311111 
0.573333 
0.835556 
1.124444 



1 .000000 
0. 907304 
0. 816445 
0.742125 
0. 670838 



COEFFICIENTS OF THE POWER SERIES EXPANSION 
Y(X1 = B(1H-B(2)*X + B(3) vX**2 . . . 



( 1) 


= 1.01512D 00 B ( 2 ) = 


-3. 186320-01 






ESTIMATES OF ERROR 


FOR THE COEFFIC I 


ENTS 


RRB ( 


1 ) = 3 . 


256D-02 ERR3( 


2) = 4. 701D-02 




SUM 


SC DEV 


= 1 . 7 55 D— 04 


F-RATIO = 


4.595D 01 


CHISQ= 1 


.0590- 03 


DEG OF FREEDOM = 


3. 


I 


X ( I ) 


F2 ( I ) 


Y ( I ) 


DELY ( I ) 


I 


X ( I ) 


F2( I ) 


Y ( I I 


DELY(I) 


1 


G.102 


1 . 00 0000 0 D 00 


9.82551 5 60— 0 1 - 


1. 7446425D-02 


2 


0.311 


9. 07303 70D-01 


9. 1599281 D-01 


8 .68911 2 0D-03 


3 


0.573 


8 .16444990-01 


8. 3 2 440 3 2D- 01 


1. 5995330D-02 


4 


0.836 


7 .42124610-01 


7 .468b783D— 01 


6 . 7 63 2205D- 0 3 


5 


1.124 


6. 70837 7 ID— 0 1 


6. 5683847D-01 - 


1 .39992370-02 
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TEMPERATURE AND RESISTIVITY VALUES 
DATA SET C 



COOLING CONDITIONS 



THE 

THE 

THE 

THE 

THE 



NUMBER OF DATA POINTS IS 21 

MAXIMUM DEGREE OF FIT IS 1 

Y-AXIS SCALE IS 0.10 IN DEGREES PER INCH 

X-AXIS SCALE IS 1.00 IN OHM-CM PER INCH 

INITIAL TEMPERATURE (TO) WAS 22.50IN DEG C 



DIMENSIONAL VALUES 



NONDIMENS IONAL VALUES 





TEMPFRATURE 


RESISTIVITY 


THETA 


RHG-PRIME 




( DEGR£ES CJ 


(OHM-CM) 








23.5 


3463.8 


0.044444 


1.000000 




25.2 


3428.4 


0.120000 


0. 989780 




25.8 


3359.8 


0.146667 


0.969975 




26.2 


3359.8 


0. 164444 


0.969975 




27.1 


3326.5 


0.204444 


0.960361 




28.4 


3261.8 


0.262222 


C .941683 




30.0 


3199. 6 


0.333333 


0.923725 




31.0 


3169.4 


0.377778 


0.915007 




31.6 


3082.1 


0.404444 


0.889803 




33.3 


3026.4 


0.480000 


0. 873723 




35.3 


2921.0 


0.566889 


0.843293 




37.2 


2822.7 


0.653333 


0.814914 




39.2 


2730.7 


0.742222 


0. 788354 




40.8 


2644.6 


0.813333 


0.763497 




46.7 


2415.9 


1.075556 


0.697471 




49 .5 


2315.8 


1.200000 


0.668572 




52.0 


2253.5 


1.311111 


0 .650586 




53.9 


2194.5 


1.395556 


0.633553 




55 .6 


2152.2 


1.471111 


0.621341 




57. 5 


2085.2 


1.555556 


0.60 1998 




61.0 


1986.3 


1.711111 


0. 573445 




COEFFICIENTS 


OF THE POWER SERI 


ES EXPANSION 






Y ( X ) = 


B ( 1 ) + B ( 2 ) <' ; X+ B ( 3 ) ‘> ,: 


X**v 2 . . . 




R ( ) ) 


= 1.00608D 00 


B { 2) =-2. 684730- 


01 






ESTIMATES OF 


ERROR FOR THE CO 


EFFICIENTS 




FRRR ( 


1 ) =2 . 339D-02 


ERR6( 2 ) = 2 . 634D- 


02 




SUM 


SO DEV = 1 .7330-04 F-RATIO = 1.039D 


02 


CHI SQ= 5 .0500-03 


DEG OF FRE 


EDOM = 19. 





1 
1 

2 
3 
A 

5 

6 

7 

8 
9 

1 0 
I 1 
12 

13 

14 
1 5 
16 
17 
1 8 

19 

20 
21 



xm 

0.044 
0.120 
0.147 
0. 164 
0.204 
0.262 
0.333 
0.378 
0.404 
0.480 
0.569 
0.653 
0.742 
0.813 
1.076 
1.200 
1.311 
1.396 
1.47 1 
1.55 6 
1.711 



F-2 ( I ) 

l.OOOOOOOD 00 
9. 8978001C-01 
9/69975 17D— 0 1 
9/69975170-01 
9. 6036145D-01 
9.416825 5D-01 
9. 23725390-01 
9. 1500664D— 01 
8.898031 ID-01 
8. 737225 OD— 01 
8^43293490-01 
8. 14914260-01 
7. 8835383D— 01 
7. 63496740-01 
6. 9747 099D-01 
6 /6 85 7 209 D- 01 
6 .50586060-01 
6. 3355275D-01 
6.21340 7 2D-01 
6.019978 ID-01 
5. 73 4453 5D-01 



Y ( I ) 

9.941484 6D-01 
9. 73 86386D-01 
9. 66704590-01 
9. 61931740-01 
9. 51 19283D-01 
9. 35681070-01 
9. 16 58968 D— 0 1 
9.04657560-01 
8. 97498290-01 
8. 7721369 D— 01 
8.53349450-01 
6. 30678420-01 
8. 068I418D-01 
7 .87722790-01 
7.1 732328 D-0 1 
6. 83913340-01 
6. 5408305 D- 01 
6. 31 4 1202 D-01 
6. 11127410-01 
5 .8845639 D-01 
5. 4669397D-01 



DEL Y( I ) 

-5 . 85 1 5376D- 03 
-i . 591 61 52D-02 
-3. 27058510-03 
-8 .04343310-03 
-9. 1686220D-03 
-6. 001470 3D— 03 
-7 . 1357053D-03 
-1. 0349077D-02 
7.69 51848D- 03 
3 .491 1861 D-03 
1. 0055955D-02 
1.57641640-02 
1 . 8460349D-02 
2. 4226050D-02 
1 .98522959-02 
1. 53412560-02 
3. 496983 6D- 03 
-2 . 1407282D— 03 
-1. 021331 10-02 
-1 .35416200-02 
-2 .67513830-02 
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TEMPERATURE AND RESISTIVITY VALUES 
DATA SET D 



HEATING CONDITIONS 



THE NUMBER OF DATA POINTS IS 5 

THE MAXIMUM DEGREE OF FIT IS 1 

THE Y-AXIS SCALF IS 0.10 IN DEGREES PER INCH 

THE X-AXIS SCALE IS 1.00 IN OHM-CM PER INCH 

THE INITIAL TEMPERATURE (TO) NAS 19.00IN DEG C 



DIMENSIONAL VALUES 
TEMPERATURE RESISTIVITY 
(DEGREES C) (OHM-CM) 



2.4.7 

32.2 

40.3 
46.5 

54 .8 



4271.7 
3551.6 
3182 .6 

2610. 8 
2496.5 



NONDIMENS IONAL VALUES 
THETA RHO-PRIME 



0.300000 

0.694737 

1.121053 

1.447368 

1.884211 



1.000000 
0. 831425 
0.745043 
0.658005 
0.584428 



COEFFICIENTS OF THE POWER SERIES EXPANSION 
Y(X) = 6(1 >+6(2 KX + b( 3 )->X**2.. . 

R( 1) = 1 . 04 268 D 00 B( 2 ) =-2 . 55 9960-0 1 

ESTIMATES OF ERROR FGR THE COEFFICIENTS 

F RRB ( 1) = 7. 1670-02 ERRB( 2 ) = 5 . 862D-02 
SUM SO DEV - 6.3490-04 F-RATIC = 1.907D 01 



CHISO= 2 


.9800-03 


DEG GF FREEDOM 


= 3. 


I 


X( 1 ) 


F2( I ) 


YU ) 


DELY(I ) 


1 


0.30 3 


t . OOOOOOOD 00 


9. 6 5 882 5 2D- 01 


-3. 41 174820-02 


2 


0.695 


8.3142 543 D-01 


8.64 8 3137 C-01 


3 .34059430-02 


3 


1.121 


7.:45G429 6D-01 


7 . 556961 4D— 01 


1. 0653179D-02 


4 


1 .447 


6 2580050 ID-01 


6. 7216052 D-01 


1. 41555140-02 


5 


1. 884 


5. S442775D— 01 


5.6033059D— 01 


-2 .40971560-02 
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TEMPER AT JRE AND 
DATA 



RESI STI VI TY VALUES 
SET D 



COOLING CONDITIONS 



THE 


NUMBER CF DATA 


POINTS IS 


19 




THE 


MAXIMUM DEGREE 


OF FIT IS 


1 




THE 


Y-AXIS SCALE IS 


0.10 IN 


DEGREES PER INCH 


THE 


X-AXIS SCALE IS 


1.00 IN 


CHM-CM PER INCH 




THE 


INITIAL TEMPERATURE (TO WAS 


19 . 00 IN DEG C 






DIMENSIONAL VALUES 


NGNDIMENSIONAL 


VALUES 




TEMPERATURE 


RESISTIVITY 


THETA 


RHO-PRIME 




(DEGREES C) 


(OHM-CMJ 








24. 1 


4166 .1 


0 .268421 


1.000000 




24.0 


4166. 1 


0.263158 


1.000000 




25.5 


4065.7 


0.342105 


0.975901 




25.7 


4017.2 


0 .352632 


0.964259 




2fl . 2 


3923.7 


0.484211 


0.941816 




29 .9 


3791.3 


0.573684 


0. 910036 




31.0 


3749. 1 


0.631579 


0.899906 




33.5 


3551.6 


0.763158 


0.852500 




35 .7 


3407.8 


0.873947 


0. 817983 




36.3 


3442.7 


1.015789 


0.826360 




42.5 


3275.4 


1.236842 


0. 786203 




46.. 4 


3066.7 


1 .442105 


0.736108 




47.6 


3066.7 


1.515789 


0.736108 




50 .4 


2498.0 


1.652632 


0. 599602 




52.0 


246 1.5 


1.736842 


0.590840 




54.7 


2408.6 


1.876947 


0. 578143 




55 .8 


2358.0 


1.936842 


0.565997 




55.3 


2189. 2 


1.910526 


0 .52p479 




57.8 


2120.3 


2. 042105 


0. 508941 




COEFF ICI ENTS 


OF THE POWER SERIES EXPANSION 






Y( X 1 = 


B ( 1 ) +5 ( 2 ) ■ X + 8 ( 


3 )*X*:*2. . . 




R ( 


1)= 1, 07 062 D 00 


B ( 2) =-2. 64046D-01 






ESTIMATES OF 


ERROR FOR THE 


COEFFICIENTS 




F RR8 { l»=5.741D-02 


FRRB( 2 ) = 4 . 532D- 02 




SUM SO DEV = 7.197 D— 04 


F-RATIO = 3.395D 


01 


CHI SQ= 2.0390-02 


DEG OF 


FREEDOM = 17. 





I 


X ( I J 


F? ( I } 


Y ( I ) 


DEL Y ( I ) 


1 


0.268 


1 . 000000 OD 00 


9. 9973970D-01 


-2. 603 00 GOD- 04 


2 


0.263 


1. OOOOOOOD 00 


1.0011294D 00 


1. 12941 360-03 


3 


0.342 


9.75900720-01 


S .8028370D-01 


4 . 38 29747D-0 3 


4 


0.353 


9 . 642591 4D- 01 


9. 7750427D— 01 


1. 3245129D-02 


5 


0.484 


9.4181609D— 01 


9. 42761 4 ID— 01 


9. 45320630^04 


6 


0.574 


9. 10035760-01 


9. 19 13626D-01 


9 . 10049 7 5 D-03 


7 


0.632 


8.9990639D— 01 


9. 0384940D-01 


3. 94301 64D- 03 


8 


0.763 


8. 52493940-01 


8.6910654D— 01 


1 . 6606603D- 02 


9 


0.879 


6. 179832 50-01 


8. 3853283D— 01 


2 .05495800-02 


1 0 


1.016 


8 226 3 6 03 9 D— 01 


8. 0240025D-01 


-2. 39601350-02 


1 1 


1.237 


7-86202920-01 


7 .44 0322 4 D-0 1 


-4. 21 7 06 79 D- 02 


12 


1.442 


7.361081 10-01 


6. 89 83338D-01 


-4.62 747 29D— 02 


13 


1 .516 


7 .361081 ID-01 


6.70377 38D-01 


-6.5 73 C 73 ID- 02 


14 


1.653 


5- 996015 5D— 01 


6. 342448 0D-01 


3 .46432590-02 


1 5 


1 . 73 7 


5». 9084035D-01 


6. 1200937D-01 


2. 1 169020D-02 


16 


1 .879 


5 278 1 42 5 3 D-0 1 


5. 74487 08 D- 01 


-3. 655 54 31 D-03 


1 7 


1.937 


5. 65996980-01 


5. 592 002 30-01 


-6 .79675000-03 


18 


1.911 


5- 2547947D-0 1 


5. 6614880D-01 


4. 06693330-02 


1 9 


2.042 


5- 0894122 D— 01 


5 .3140594 D— 0 1 


2 . 246 47 2 ID-02 
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TEMPERATURE AND RESISTIVITY VALUES 
DATA SET E 



HEATING CONDITIONS 



THE NUMBER OF DATA POINTS IS 5 

THE MAXIMUM DEGREE OF FIT IS 1 

THF Y-AXIS SCALE IS 0.10 IN DEGREES PER INCH 

THE X-AXIS SCALE IS 1.00 IN OHM-CM PER INCH 

THE INITIAL TEMPERATURE (TO) WAS 22.50IN DEG C 



DIMENSIONAL VALUES 
TEMPERATURE RESISTIVITY 
(DEGREES C) ( OHM-CM J 



25.7 

35.7 
41 .0 

44.7 
49.5 



3451.0 

3105.5 
2894.7 

2710.6 

2529.6 



NGNDIMENSIONAL VALUES 
THETA RHO-PRIME 



0.142222 
0.506667 
0 .822222 
0.906667 
1.200000 



1 .COOOOO 

0. 899804 
0.838800 
0.785453 
0. 733005 



COEFFICIENTS OF THE POWER SERIES EXPANSION 
Y ( X ) = B ( 1 I + B ( 2 ) 4 X + B ( 3 )4'X**2. • • 

R ( 1) = 1.04182D 00 B ( '2) = -2. 54679D- 0 1 

ESTIMATES OF ERROR FOR THE COEFFICIENTS 

F R P, B ( 11 = 1 .6920-02 cRRB( 2) = 2.036D-02 



SUM 


SO DEV 


= 3. 281D-05 


F-RATIC 


= 1.5650 02 


CHI SG= I 


.8930-04 


DEG OF FREEDOM 


= 3. 


I 


X( I ) 


E2< I) 


Y ( I ) 


DELY ( I ) 


1 


0.142 


UOOOOOOOD 00 


1. 0055945D 00 


5. 59454790-03 


2 


0.587 


8.9988409 D-G 1 


8. 92403670-01 


-7.4804221D-03 


3 


0.822 


8. 3880035D— 01 


8. 3241250D— 01 


-6 . 3 8 7843 9 D— 03 


4 


0.987 


7.8545349D-01 


7. 9053188 D- 01 


5. 07838 71D-03 


5 


1.200 


7.3300493 D— 01 


7 .3620026D— 01 


3 . 19533100-03 
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TEMPFRATURE AND RESISTIVITY VALUES 
DATA SET E 



COOLING CONDITIONS 



THE NUMBER OF DATA POINTS IS 9. 

THE MAXIMUM DEGREE GF FIT IS 1 

THE Y-AXIS SCALE IS 0.10 IN DEGREES PER INCH 

THE X-AXIS SCALE IS 1.00 IN OHM-CM PER INCH 

THE INITIAL TEMPERATURE (TO) WAS 22.50IN DEG C 



DIME NS I 
TEMPERATURE 
(DEGREES C) 

24.3 
26. .4 

27.5 

29.5 

30.5 

32.0 

33.0 

33.3 

63.5 



ONAL VALUES 
RESISTIVITY 
(OHM-CM ) 

3796.6 

3674.0 

3559. 0 
3522.3 
3416 .4 

3382.6 

3316.8 

3284 .9 

2627. 1 



NCNDIMENS I ONAL VALUES 
THETA RHG-PRI ME 



0.080000 

0.173333 

0.222222 

0.311111 

0.355556 

0.422222 

0.466667 

0.480000 

1.822222 



1.000000 
0.967708 
0 .937418 
0. 927751 
0.899858 
0.890955 
0. 873624 
0.865222 
0.691961 



COE FF I C I.E NT S OF THE POWER SERIES EXPANSION 
Y ( X ) = B ( 1 ) +B ( 2 ) i X*- B ( 3 ) - X**2. . . 



R ( 1)= 9. 73 366D- 01 B( 2) =-l. 62 877D-01 

ESTIMATES OF ERROR FOR THE COEFFICIENTS 



RRB ( 


11 = 3. 


445D-02 ERR3( 


2) =5. 010D-02 




SUM 


SO DEV 


= 4.706D-04 


F-RAT 10 


= 1.057D 01 


CH I SQ= 4 


.65.78-03 


DEG OF FREEDOM 


= 7 . 


I 


X( I) 


F2 ( I ) 


Y ( I ) 


DEL Y( I ) 


1 


0. 080 


1 . OOOOOOOD 00 


9.6033601 D— 01 


-3.96639890-02 


2 


0.173 


9. 677 0795 D- 01 


9.451341 60—01 


-2. 25737688-02 


3 


0.222 


9.3741 76 9D-01 


9.37171320-01 


-2 .46371190-04 


4 


0.311 


9.277511 5 D-0 1 


9. 226933 60-01 


-5 .05776030-03 


5 


0.356 


8 . 998 57 77D-0 1 


9 . 1 545442u- 0 1 


1. 55966520-02 


6 


0.422 


8 .'906 55 07 D— 01 


9.0459597 D— 0 1 


1 .364090 40— 0 2 


7 


0.467 


8. 73623 7 7D- 01 


8 . S 7357 0 OD— 01 


2.3733234D-02 


8 


0.480 


8 .65221510-01 


8.951853 ID- 01 


2 . 9963799D- 02 


9 


1. 822 


6.:51 961 2 3D-01 


6. 76568530-01 


-1 .53927010-02 
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TEMPERATURE AND RESISTIVITY VALUES 
DATA SET F 



HEATING CONDITIONS 



THE NUMBER OF DATA POINTS IS 5 

THE MAXIMUM DEGREE OF FIT IS I 

THE Y-AXIS SCALE IS O.IO IN DEGREES PER INCH 

THE X-AXIS SCALE IS 1.00 IN OHM-CM PER INCH 

THE INITIAL TEMPERATURE (TO) WAS 22.50IN DEG C 



DIMENS I 
TEMPERATURE 
(DEGREES C) 

24.2 

30.0 
33.. 0 
44. 8 

49.0 



ONAL VALUES 
RESISTIVITY 
(OHM-CM) 

3298.8 
3009. 1 

2559.3 

2381.3 

2197.8 



NGNDIMENSIONAL VALUES 
THETA RHO-PRIME 



0.075556 
0. 333333 
0 .686389 
0.991111 
1.177776 



1 .000000 
0.912180 
0. 775828 
0 .721669 
0.666242 



COEFFICIENTS OF THE POWER SERIES EXPANSIGN 
Y ( X ) = B ( 1 )+B(2)*X + 6(3)*X4*2.. . 

R( 1)= 1.011910 00 B ( 2)--3.01045D— 01 

ESTIMATES OF ERROR FOR THE COEFFICIENTS 



F R R B ( I ) = 3 . 607 D-0 2 E R R S ( 
SUM SO DEV = 2. 179D-04 



CHI S0= 1 


.3600-03 


I 


X( I) 


F2 ( I ) 


1 


0.076 


1.00000000 00 


2 


0.333 


9. 12180190-01 


3 


0.689 


7.75827570-01 


4 


0.991 


7.21868560-01 


5 


1.178 


6 .66242270-01 



2 ) = 4. 687D-02 

F-RAT ID 
DEG OF FREEDOM 

Y ( I) 

9. 89161 02D-C1 
9.1 155822D— 01 
6.0451988 D— 0 1 
7. 13 5373 CD- 01 
6. 5734217D-01 



= 4. 1250 01 
= 3. 

DEL Y (I ) 

-1 . 083 8983D-02 
-6. 21964870-04 
2.86923100-02 
-8 . 331 26 16 D-0 3 
-8.90010070-03 
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TEMPERATURE AND 
DATA 



RESISTIVITY 
SET F 



VALUES 



COOLING CONDITIONS 



THE NUMBER OF DATA POINTS IS 12 

THE MAXIMUM DEGREE OF FIT IS 1 

THE Y-AXIS SCALE IS O.IO IN DEGREES PER INCH 

THE X-AXIS SCALF IS 1.00 IN OHM-CM PER INCH 

THE INITIAL TEMPERATURE (TO) WAS 22.50IN DEG C 



DIMENSIONAL VALUES 


NONDIMENS IONAL 


VALUES 


TEMPERATURE 


RES1 STI VI TY 


THETA 


RHO-PRIME 


(DEGREES C) 


(OHM-CM ) 






24.0 


3630.9 


0.066667 


1 .000000 


25.0 


3581.6 


0.111111 


0.986422 


25.5 


3574.1 


0.133333 


0.984356 


27.9 


3465. 7 


0.240000 


0.954502 


29.8 


3330.9 


0.324444 


0.917376 


33.5 


3206.2 


0.488889 


0.883032 


38.2 


2906.9 


0.697778 


0.800600 


43-. 0 


2743.9 


0.911111 


0. 755708 


51.0 


2449.4 


1.266667 


0.674599 


5 5..0 


2348.6 


1.444444 


0.646837 


60.5 


2183 .8 


1 .688889 


0.601449 


65.5 


2077.7 


1.911111 


0.572227 


COEFFICIENTS 


OF THE POWER S 


ERIES EXPANSION 




Y( X ) = 


B ( 1 ) +B ( 2)*X+B(3)*X*32... 




R< L)= 1. 00 2 ODD 00 


B( 2) =-2.420020-01 




ESTIMATES OF 


ERROR FOR THE 


COEFFICIENTS 




F RR3 ( 1 ) = 3. 17-5D-02 


ERRB( 2 ) =3. 181D— 02 




SUM SO DEV = 3 .350 0- 04 


F-RATIO = 5.786D 


01 


CHI SG= 5 .6070—03 


DEG OF 


FREEDOM = 10. 





I 


X ( I ) 


F? ( I ) 


Y ( I) 


DELYCI ) 


1 


0.067 


l . OOOOOOOD 00 


9 .85 86368 D-01 


-1 .41363223-02 


2 


0.111 


9 .- 8642210 D— 0 1 


9.7 51030 ID— 01 


-1. 1314087D— 02 


3 


0.133 


9 .1843565 0 D— 0 1 


9.6973018 D— 0 1 


— 1 .462 63 1 oQ— 02 


4 


0.240 


9- 54501 64D— 01 


. 9. 43 9 1 65 8 D— 01 


-1 . 0585056D— 02 


5 


0.324 


9 . 173758 6D-0 1 


9. 23480 82 D- 01 


6. 10496100-03 


6 


0.489 


8.8-303176 D— 01 


8.8368486 D— 01 


6 .53100990-04 


7 


0.698 


8.1006004 0D-01 


8. 33133230-01 


3.25328270-02 


8 


0.911 


7.I5570795D-01 


7.8 1 5060 3D— 0 1 


2. 57980840-02 


9 


1.267 


6.745985 8 D-01 


6 .954607 ID-01 


2 .08621270-02 


10 


1.444 


6-. 4683687D- 01 


6. 5243805D-01 


5.601 17 74D— 03 


1 1 


1.689 


6*0144863 D— 0 1 


5. 9328189 D-01 


-8. 166 7 866D- 03 


12 


1.911 


5. 722272 7D-01 


5. 39 503560-01 


-3 .2723708 D-02 
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TEMPFRATURE AND RESISTIVITY VALUES 
DATA SET G 



HEATING CONDITIONS 



THE NUMBER OF DATA POINTS IS 5. 

THE MAXIMUM DEGREE OF FIT IS 1 

THE Y— ,AX I S SCALE IS 0.10 IN DEGREES PER INCH 

THE X-AXIS SC A L c IS 1.00 IN OHM-CM PER INCH 

THE INITIAL TEMPERATURE (TO) WAS 15.00IN DEG C 



DIMENSI 
T EMPERATURE 
(DEGREES C) 

25 .0 

30.3 
39.0 
47 .3 

51.4 



ONAL VALUES 
RESISTIVITY 
(OHM-CM) 

3134.0 

2022.8 

2439. 1 
2188.5 
2031. 9 



NONDI MENSI ONAL VALUES 
THETA RHO-PRIME 



0.666667 

1.020000 

1.600000 

2.153333 

2.426667 



1.000000 

0.900702 

0.778271 

0.698309 

0.648341 



COEFFICIENTS OF THE POWER SERIES EXPANSION 
Y( X) = 3( 1 )+B( 2)4X+B( 3)--X?*2. . . 



B( 1)= 1.11070D 00 B ( 2) =-l .942230-01 

ESTIMATES OF ERROR FOR THE COEFFICIENTS 



FRRB ( l ) = 4. 886D- 02 ERR6( 2)=2.862D-02 
SUM SO DEV = 2.156D-04 F-RATIO = 4. 6050 01 

CHI SQ= 1.27 oD-03 DEG OF FREEDOM = 3. 



1 

2 

3 

4 

5 



0 .667 
1.020 
1.600 
2.153 
2.427 



1 .:0 0 0 0 0 0 0 D 00 
9* 00 701 98 D— 01 
7. .782705 80- 0 1 
6.9330887 D— 0 1 
6. 483407 3D— 01 



9. 8 122036 D- 01 
9. 125947 4 D— 01 
7. 99945150-01 
6.92 4748 5D— 01 
6.3938711 D— 0 1 



-1 . 8 77 9641D-02 
1 . 189 27 66D-0 2 
2. 1 674569D— 02 
-5 . 8340231D— 03 
-8 .95367 10D-03 
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TEMPERATURE AND RESISTIVITY VALUES 
DATA SET G 

COOLING CONDITIONS 



THE 


NUMBER OF DATA 


POINTS IS 


10 




THE 


MAXIMUM DEGREE 


OF FIT IS 


1 




THE 


Y-AXIS SCALE IS 


0 . I 0 IN 


DEGREES PER INCH 


THE 


X-AXIS SCALE IS 


1.00 IN 


OHM-CM PER INCH 




THE 


INITIAL TEMPERATURE (TO) WAS 


15. 00 IN DEG C 






DIMENSIONAL VALUES 


NON DIMENSIONAL 


VALUES 




TEMPERATURE 
(DEGREES C) 


RESISTIVITY 

(OHM-CM) 


THETA 


RHO-PRIME 




23.0 


3416.5 


0 .533333 


1.000000 




24.7 


3253.6 


0 • 6466 6 7 


0.952320 




28 .0 


3105.5 


0.866667 


0.908971 




30.7 


2970.3 


1.046667 


0.869399 




35.5 


2732.3 


1.366667 


0. 799737 




3 9 .0 


2587.2 


1 .600000 


0.757266 




45 . 5 


2371.3 


2.033333 


0 .694073 




50.8 


2217.0 


2.386667 


0. 648910 




56.3 


2069.0 


2.753333 


0.605591 




63.2 


1917. 5 


3.213333 


0.561247 



COEFFICIENTS OF THE POWER SERIES EXPANSION 
Y(X1=3( 1)+B( 2)*X+B(3)^X=?*2... 

R( 1) = 1.04624D 00 B( 2) =-l. 620310-01 



F RRB ( 


ESTIMATES OF ERROR 
1 ) = 5. 445D- 02 ERRB( 


FOR THE COEFFICIENTS 
2 ) =2. 922D-02 


SUM 


SO DEV 


= 5 .2 360-04 


F- RAT 10 


= 3.074D 01 


CHI S0= 7 


.1220-03 


DEG OF FREEDOM 


= 8. 


I 


X( I) 


F? ( I) 


Y ( I) 


DEL Y ( I ) 


l 


0. 533 


1. 00 00 GOOD 00 


9.5982U6D-01 


-4. 0178821 D— 02 


2 


0.647 


9-. 523 1 96 3D-01 


9.41457710-01 


-1. 086 19 10D— 02 


3 


0.867 


9. 08971170-01 


9 . 0 5 8 109 9 D-0 1 


-3 . 1601795D-C3 


4 


1.047 


8 2693 985 ID— 01 


8. 76645490-01 


7 .24698020-03 


5 


1.367 


7.99 73657D-01 


8.24795710-01 


2. 5059133D-02 


6 


1.600 


7. 572662 1 D-Oi 


7.86588570-01 


2 . 972 2365D— 0 2 


7 


2.033 


6. 94072880—01 


7. 167753 3D— 01 


2 .27024450-02 


8 


2.387 


6 t 48 90 97 OD— 01 


6 . 59 5245 3D— 01 


1 * 06 14322D— 02 


9 


2. 753 


6. 0559052 D-Oi 


6.0011332 C— 0 1 


-5 .47719970-03 


10 


3.213 


5261246890—01 


5. 2557926D-01 


-3. 5667635D-02 
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TEMPERATURE AND RESISTIVITY VALUES 
DATA SET H 



HEATING CONDITIONS 



THE NUMBER OF DATA POINTS IS 5 

THE MAXIMUM DEGREE OF FIT IS 1 

THE Y-AXIS SCALE IS 0.10 IN DEGREES PER INCH 

THE X-AXIS SCALE IS 1.00 IN OHM-CM PER INCH 

THE INITIAL TEMPERATURE (TO) WAS 7.50IN DEG C 



DIMENS I 
TEMPERATURE 
(DEGREES C) 

24 .7 
30.0 

40.7 
45.4 

50.7 



ONAL VALUES 
RESISTIVITY 
( OHM-CM) 

3062.9 

2744.0 
2316.8 
2142.7 

1970.0 



NONDI MEN S I ONAL VALUES 
THETA RHO-PRIME 



2.293333 
3.000000 
4 . 42666 7 
5 .053333 
5.760000 



1.000000 
0.095883 
0. 756407 
0.699566 
0.643181 



COEFFICIENTS OF THE POWER SERIES EXPANSION 
Y ( X ) = B (1) +3 ( 2 ) ••- X + 3 ( 3 ) -> X**2 . . . 



R I 1)= 1.215940 00 B ( 2) =-l . 0 1 526D-01 

ESTIMATES OF ERROR FOR THE COEFFICIENTS 



FRRB( 1 ) = 5. 41 3D- 02 E R R B ( 21=1.2580-02 

SUM SO DEV = 1 . 566 D— 04 F-RATIO = 6.5120 01 

CH I SQ= 9.317D-04 DEG OF FREEDOM = 3. 



I X ( I ) 

1 2.293 

2 3.000 

3 4.427 

4 5.053 

5 5.760 



F? ( I ) 

1 . 000000 OD 00 
3 ..95882990-01 
7 -.56407330-01 
6.99565770-01 
6 .43181300—01 



Y { I ) 

9.3310763 D— 0 1 
9. 11362720-01 
7.6o 5 192 ID— 0 1 
7. 02896360-01 
6. 31151450-01 



DEL Y ( I ) 

-1 .68923650-02 
1.54797330-02 
1 .01118870-02 
3.33059360-03 
- 1. 20298490-02 
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TEMPERATURE AND RESISTIVITY VALUES 
DATA SET H 



COOLING CONDITIONS 



THE NUMBER OF DATA POINTS IS 
THE MAXIMUM DEGREE OF FIT IS 
THE Y-AXIS SCALE IS 0.10 IN 

THE X-AXIS SCALE IS 1.00 IN 

THE INITIAL TEMPERATURE (TO) WAS 



10 

1 

DEGREES PER INCH 
OHM- CM PER INCH 
7 .50 IN DEG C 



DIMENSIONAL VALUES 


NONDIMENS IONAL 


VALUES 


TEMPERATURE 


RESISTIVITY 


THETA 


RHO-PRIME 


(DEGREES C) 


(OHM-CM ) 






24.0 


3267.4 


2.200000 


1 .000000 


24.5 


3267.4 


2.266667 


1.000000 


, 25.0 


3267.4 


2 .333333 


1.000000 


26.5 


3147.3 


2.533333 


0.963243 


31.4 


2901.0 


3.186667 


0. 887862 


37.2 


2658.7 


3.960000 


0.813705 


44.5 


2364. 3 


4.933333 


0.723756 


49 .5 


2226.4 


5.600000 


0.681398 


65.5 


1842.6> 


7.733333 


0.563935 


70.5 


1757.4 


8.400000 


0. 537859 


COEFFICIENTS 


OF THE POWER 


SERIES EXPANSION 




Y ( X ) = 


B ( 1 )+3(2)*X+B(3)vX**2... 




R ( 1»= 1. 15284D 00 


B ( 2)=-7.77955D-02 




ESTIMATES OF 


ERROR FOR TH 


E COEFFICIENTS 




FRRB ( 1 ) = 7.0 09D-02 


ERRB( 21=1. 450D— 02 




SUM SO DEV = 7 . 9990-04 


F-RATIG = 2.880D 


01 


CHI SQ= 1.CS9D-02 


DEG OF 


FREEDOM * 8. 





I X ( I ) F? ( I ) 

1 2.200 1 . OOOOOOOD 00 

2 2.267 1.00000000 00 

3 2.333 1*00000000 00 

4 2.533 9. 6324295D-0 1 

5 3.187 6*673619 1Q-G1 

6 3.960 8. 1370509D— 01 

7 4.933 7. 237553 9D- 0 1 

8 5.600 6.'8139805D— 01 

9 7.733 5j. 6393463D-01 

10 8.400 5 . 378588 5D— 01 



Y ( I ) 

9. 81687270-01 
9.7to50090D— 01 
9.7131453 D-01 
9. 55 7 5544D— 01 
9. 049 29 05 D-01 
8. 447672 OD— 01 
7. 6904626D-01 
7. 1718260 C— 01 
5. 5121839D-01 
4. 9935523D- 01 



DEL Y ( 1 ) 

•1 . 831 2735D— 02 
•2.34991010-02 
-2 .86854670 -02 
-7.487 51 03 D- 03 
1 . 70671400-02 
3 . 10621160-02 
4.^52903680-02 
3 . 5784546D-02 
-1 .27157390-02 
-3. 85036200-02 
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SAMPLE DATA 
SPHERICAL PROBE 

DATA SET I 



T 


I 


TAU 


R 


R/RO 


R R/RO 


R R/RO 


( SFC ) 


(AMPS) 




(CM) 




(CM) 


( CM) 








R— 27 




ROCHROME 


R-37 


30. 0 


C.050 


0.78 


1 . 05 


1.62 






60 .0 


0.050 


1.56 


1.20 


1.85 


NOT 


NOT 


120.0 


0.050 


3.12 


1.45 


2.23 






240.0 


0.050 


6.25 


1.85 


2.85 


USED 


USED 


480 .0 


0.050 


12-.50 


2.00 


3.08 






960.0 


0.050 


25.00 


2.30 


3.54 







THETA/GO 0.1260 

RESISTIVITY (OHM-CMJ 2420.0 

ISOTHERM (DEGREE Cl 27.9 

THE INITIAL TEMPERATURE WAS 19.2 DEGREES C. 

THF PROBE RADIUS (ACTUAL) WAS 0.241 CM. 

THE PRGBE RADIUS IN THE PICTURES WAS 0.650 CM. 



SPHERICAL PROBE 



DATA SET J 



T 


I 


TAU 


R R/RO 


R 


R/RO 


R R/RO 


(SFC ) 


( AMPS ) 




(CM) 


(CM) 




(CM) 








R— 27 


ROCHROME 


R-37 


30.0 


0.06 0 


0.78 


1.10 2. SO 


0. 90 


2.37 


0.0 0.0 


60.0 


0.060 


1.56 


1.40 3.68 


1.10 


2.90 


0.50 2.37 


120. 0 


0. 060 


3.12 


1.65 4.34 


1.40 


3.68 


1.00 2.63 


240.0 


0.060 


6.25 


2.15 5.66 


1.70 


4.47 


1.25 3.29 


480.0 


0. 060 


12.5 


2.65 6.97 


2.05 


5.40 


1.45 3.82 


960.0 


0.060 


25.00 


3.30 8.68 


2. 55 


6. 71 


1.70 4.47 


THF TA 


/GO 




0.0200 


0.0282 


0.0572 


RESISTIVITY (OHM-CM) 


3605.0 


3520 


. 0 


3190.0 


ISOTHERM (DEGREE C) 


27.9 


29. 


9 


37.0 



THF INITIAL TEMPERATURE WAS 23.0 DEGREES C. 
THE PROBE RADIUS. (ACTUAL) WAS 0.241 CM. 

THF PROBE RADIUS IN THE PICTURES WAS 0.38 CM. 
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SAMPLE DATA 

CYLINDRICAL PROBE (L/D=4) 



TAU 



DATA SET K 
R/RO 



R/RO 



RESISTIVITY (GHM-CM) 
ISOTHERM (DEGREE C) 



3465 .7 
27.9 

THF INITIAL TEMPERATURE WAS 22. 

THE PROBE RADIUS (ACTUAL) WAS O.I2I CM. 

THE PROBE RADIUS IN THE PICTURES WAS 0.19 CM. 
THF PROBE LENGTH WAS 0.96 CM. 



3324.9 
29.9 
DEGREES C. 



R/RO 



(SFC) 


(AMPS) 




(CM ) 




(CM) 




(CM) 








R- 


27 


ROCHROME 


R-37 


30.0 


0.030 


3.11 


0.70 


3 .68 


0.0 


0.0 


NOT 


60.0 


0.030 


6.23 


0.90 


4. 74 


0. 50 


2.63 




120 .0 


0.030 


12.45 


1 .15 


6.05 


0. 80 


4.21 


USED 


240. 0 


0. 030 


24.91 


1.42 


7 .47 


1.10 


5 .79 




THFTA/GG 




0 


.367 


0 . 


501 





CYLINDRICAL PROBE (L/D=10) 
DATA SET L 



T 


I 


TAU 


R 


R/RO 


R 


R/RO 


R R/RO 


(SFC) 


(AMPS 1 




(CM) 




(CM ) 




(CM) 








R-27 




ROCHROME 


R-37 


60. 0 


0. 090 


6.23 


0.80 


4.00 


0 .50 


2.50 


NOT 


120.0 


0.090 


12.45 


1.20 


6. 00 


0. 90 


4.50 


USED 


240 .0 


0.090 


24 .91 


1 .55 


7.75 


1.30 


6.50 





THF T A / GO 0.-357 0.452 

RESISTIVITY i OHM -CM ) 3465.7 3324.9 

ISOTHERM (DEGREE 0) 27.9 29.9 

THF INITIAL TEMPERATURE WAS 20.4 DEGREES C. 
THF PROBE RADIUS (ACTJaL) WAS 0.121 CM. 

THF PROBE RADIUS IN THE PICTURES WAS 0.200 CM 
THF PROBE LENGTH WAS 2.42 CM. 



CYLINDRICAL PROBE (L/D=20) 
DATA SET M 



T 


I 


TAU 


R 


R/RO 


R 


R/RQ 


R 


R/RO 


(SFC) 


(AMPS) 




(CM) 




(CM) 




(CM) 










R- 


27 


ROCHROME 


R- 


37 


30.0 


0. 150 


3.11 


0.70 


4.31 


0.48 


2. 95 


C. 0 


0. 0 


60.0 


0.150 


6 .23 


1 .20 


7 .39 


0.80 


4.92 


0.0 


0.0 


120.0 


0. 150 


12.45 


1.62 


9. 97 


1.15 


7.08 


0.45 


2.77 


240.0 


0.150 


24.91 


2.20 


13.54 


1.62 


9.97 


0.70 


4.31 


THFTA/GO 




0 


.356 


0.489 


NOT 


RESISTIVITY (OHM-CM) 


3551.7 


3480 


.0 


USED 


ISOTHFRM (DEGREE C) 


27.9 


29 


.9 







THF INITIAL TEMPERATURE WAS 22.5 DEGREES C. 
THF PRO B F RADIUS (ACTUAL) WAS 0.121 CM. 

THF PROBE RADIUS IN THF PICTURES WAS 0.163 CM 
THF PRCBE LENGTH WAS 4.84 CM. 
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